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Résumé 


La construction de la nouvelle laverie de Benguérir entre dans la perception d’augmenter la 
production des phosphates afin de répondre aux attentes des marchés internationaux, la 


laverie en phase d’exécution est constituée de : 


vV Station de lavage 
V Station de flottation 
v  Décanteur 


v Digue d'épandage 


Le projet étudié au sein du groupe OCP consiste à étudier l’écoulement de la boue à 
l’intérieur de la digue d’épandage afin d’optimiser le taux de stockage de solide et le taux de 


récupération d’eau. 


Pour ce faire, une caractérisation du produit boue est nécessaire, afin de modéliser son 
comportement à l’intérieur de la digue et évaluer le taux de récupération d’eau pour plusieurs 


scenarios (agir sur la pente du bassin, introduire des chicanes...). 


Le but est de trouver la configuration optimale du bassin en terme de stockage et de 
récupération d’eau et sortir avec un planning d’utilisation de la digue sans oublier de 


dimensionner les ouvrages constituants le bassin. 
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ABSTRACT 


ECOLE MOHAMMADIA D'INGÉNIEURS 


The construction of Benguerir's new beneficiation plant is part of the perception of increasing 
the production of phosphates in order to meet the expectations of the international markets. 


The beneficiation plant at the execution phase consists of: 


v Wash station 

v Flotation station 

v Thickener 

vV Spreading dyke (tailings storage facility) 


The project studied within the OCP group consists of studying the flow of sludge inside the 
spreading dike (tailings storage facility) to optimize the solid storage rate and the water 


recovery rate. 


To do this, a characterization of the sludge product is necessary, in order to model its behavior 
inside the dike and evaluate the water recovery rate for several scenarios (act on the slope of 


the basin, introduce baffles ...). 


The goal is to find the optimal configuration of the pond in terms of water storage and 
recovery and to come out with a schedule of use of the dike (tailings storage management) 


without forgetting to size the structures constituting the tailings storage facility. 
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Introduction : 


La qualité de la roche extraite dans les différents sites du Groupe lui attribue une place 
privilégiée sur l’échiquier mondial des engrais. La concentration de P205 dans la roche de 
phosphate produite au Maroc varie entre 5% et 45%. La concentration moyenne s’élève à 
31,5%. Selon les estimations du Groupe, les réserves de Khouribga représentent environ 44% 
du total des réserves disponibles au Maroc, alors que les réserves des sites de Gantour et 
Boucraû représentent respectivement 36% et 2%. 

Pour son programme de développement industriel 2008-2025 dont l’ambition est de 
servir 50% de la demande incrémentale du marché des engrais, le Groupe OCP prévoit de 
doubler la production de phosphate pour atteindre 57,4 millions de tonnes par an contre 28 
millions de tonnes actuellement, d’où la construction de la nouvelle laverie de Benguerir 
couvrant ainsi une partie du gisement Gantour de capacité 16 Milliards m° en matière de 
phosphate. 

Le but derrière cette étude est de trouver la configuration optimale du bassin de 
rétention de cette nouvelle laverie, en terme de stockage de boue et de récupération d'eau et de 
sortir avec un planning d'utilisation sur la durée de vie du bassin ainsi que le 
dimensionnement des ouvrages constituants la digue d'épandage. 

Pour ce faire, le premier chapitre concerne une présentation globale des processus de 
traitement de phosphates dans cette laverie, le deuxième introduit la digue d'épandage, le 
troisième et le quatrième chapitre contiennent une étude théorique du comportement de la 
boue à l'intérieur du bassin et une simulation de son écoulement en utilisant un logiciel appelé 
RIFT TD pour plusieurs scenarios de configuration de la digue afin de trouver le scenario 
optimal et le planning d'utilisation correspondant. Le dernier chapitre entame le 
dimensionnement des composantes du bassin, à savoir la stabilité du corps de la digue en 


utilisant geoslope et la conception des chicanes sur CYPE. 
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I. Gisements phosphatés marocains : 


Le Maroc est caractérisé par quatre gisements phosphatés, ces derniers sont déposés 
au niveau des golfs à partir du crétacé supérieur jusqu'à l’éocène, suite à une transgression 
très importante. Le phosphate déposé est de type sédimentaire sous-marin et présente 75% des 


réserves mondiales. Ces gisements sont : 


Æ Bassin D’Oulad Abdoun : 37,35 Milliards de m°. 
% Bassin des Gantour: 15,95 Milliards de m°. 

% Bassin de Meskala: 1,11 Milliards de m°. 

% Bassin Oued Eddahab : 31,35 Milliards de m°. 


Tanger 


- Bassin des Ouled Abdoun 
= Hassan des tGanmour 

3- Bassin des Meskala 
- Bassin dés Oued Eddahab 





Figure 1: Gisements phosphatés au Maroc 


Notre zone d’étude se situe dans le plateau de phosphate des Gantour, orienté Est-Ouest, de 
longueur 120 Km et de largeur de 20Km à 30km, limitée au Nord par le massif paléozoïque 
de Rehamna, au Sud par le massif ancien des Jbilet , A Est par oued Tassaout, (affluent de 


l’oued Oum-Er-rbia), à l’ouest, les collines jurassiques de mouissate. 
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Figure 2: Gisement de GANTOUR 


IT. Représentation des étapes de traitement des phosphates : 





On coupe les inférieures à 40um 


Alimentation í d 


en 
phosphates 
à partir des 





Flottation Décantation Digue d'épandage 





Grains Grains à 
de 40 à de 40 à diamètre 


200um 200um <200um 






Grains à 
Diamètre 
<10mm 


zones 
VA . 
d'extraction 





Afin d’avoir une idée sur la richesse en phosphates selon la granulométrie, on prend des 
échantillons à partir des couches de phosphates, on représente la teneur en BPL en fonction de 


la taille des grains pour savoir à quel diamètre on doit couper dans la laverie. 
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Une représentation d’un échantillon de la première couche : 


% BPL (teneur en phosphate) 
Eu 


72 





PHOSPHATE 
MARCHAND 






RS r 


40 um 180 pm 2400 um 10000 pm 





a 


Élimination 
des 





Figure 3: la teneur en BPL (Bone Phosphate of lime) en fonction du diamètre des grains 


+ Lavage: 


Pour effectuer ce processus, les grains de diamètre inférieur à 10mm reçues à partir des 


couches de phosphates en zones d’extraction sont regroupés par taille : 


Les particules appartenant à [3150um ; 10mm] : sont très pauvres en teneur en BPL, c’est 
pourquoi elles sont mises à terril et stockées sous formes de cavaliers, qui vont être utilisés 


comme matériaux de construction de la digue d’épandage. 


Les particules à diamètre inférieur à 40um : sont également pauvres en teneur en BPL, et sont 


dirigées directement vers le décanteur. 


Les particules appartenant à [200um ; 31 50m] : sont riches en BPL (environ 68%) et 
considérés comme phosphates pures. Après leur lavage, elles sont dirigées directement vers le 


stock. 


Les particules appartenant à [40um ; 200um] : sont moyennes en teneur en BPL (environ 


45%) et sont dirigées vers le block de flottation pour leur enrichissement en phosphate. 





Mémoire du projet de fin d'études À 
Æ Flottation : 
Ce block reçoit les particules en moyenne teneur en BPL, les essais ont montrées qu’ils 


appartiennent à la gamme [40um ; 200um], ce processus a pour but d’enlever les Si02 et les 


MeN et CO: afin d’avoir des phosphates pures, ces derniers sont dirigés vers le stock. 
Les matières stériles sont dirigées vers le décanteur. 


4 Décantation : 


Ce block reçoit les stériles inférieurs à 200um avec un total de solide égal à 7,42%. Après leur 
épaississement, on récupère à peu près 80% de l’eau, stocké dans un bassin qui alimente 
l’usine de production. A la sortie du décanteur, on obtient une boue avec un total de solide 


égal à 35%, stocké dans la digue d’épandage. 





Figure 4: épaississeur de boues 


+ Digue d’épandage : 


La digue d’épandage est un bassin rectangulaire étendu (1250x624 m°) avec une hauteur de 
13m, les limites du bassin sont sous forme de digue symétrique de fruits 2.6H/1V (annexe 


1.c). Le concept de récupération d’eau est le suivant : 


On achemine la boue du décanteur vers la digue, la boue s’écoule le long du bassin en 
longueur et elle se décante, après un certain moment on obtient une séparation des phases 


solide- liquide. 


Dans l’autre côté du bassin, on installe des clarinettes conduites métalliques qui collectent 
l’eau flottante et la dirige à travers une tranchée vers le bassin de reprise d’eau. (Voir plans du 


bassin en annexe 1.c). 
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Chapitre II : Présentation de la digue d'épandage de la 
nouvelle laverie de Benguerir 
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Afin de pouvoir stocker la boue provenant du circuit du lavage de phosphates et en 
extraire l’eau résiduelle, le projet de la future laverie de BENGUERIR prévoit l’aménagement 
d’une digue de rétention, dont emplacement a fait l’objet d’étude de comparaison technico- 
économique. Concernant le contexte géographique du projet, la ville de Benguerir se situe à 
environ 70 Km au Nord de Marrakech et à 57 Km à l’Ouest de la ville d’El Kelaâ des 
Sraghna. 


I. Les variantes du site de la digue : 


Fr 
Ao 





Figure 5: variantes de l'emplacement de la digue 


Variante 1 : Afin de pouvoir stocker la boue provenant du circuit de la laverie et en extraire 
l’eau résiduelle, JESA a prévu initialement, aménagement d’une digue de rétention classique 
(variante 1 sur photo), dont emplacement proposé est au Sud-Ouest de la laverie (rectangles 


jaunes), sur un terrain vierge acquis par OCP. 


Variante 2 : Succession de digues aménagées au droit des tranchées, creusées lors de 
lexploitation minière, moyennant la fermeture des cavaliers par un remblai (rectangle vert). 


L’étendu de ces digue « en parallèle » est tributaire du volume total de la boue à traiter. 


-= Ce fut l'implantation au niveau de la zone des cavaliers (variante 2) que retiennent le 
concepteur (JESA) et le Client (OCP) pour le stockage des boues issues du lavage. Ce choix 
traduit leur volonté conjointe d’apporter de optimisation en termes d’occupation des 


terrains alloués au stockage de la boue, combinée à la nécessité et l’obligation de 
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réhabiliter les anciennes mines exploitées et réduire ainsi impact de l’activité minière 


sur l’environnement. 


II. Etude du terrain d'implantation : 


l. 


Topographie : 


La topographie de la zone est relativement tabulaire avec une pente quasi-uniforme et 


régulière qui décroit du Nord vers le Sud, généralement inférieure à 10°. Dans quelques zones 


réduites, la pente peut dépasser 10°. Le nivellement de la zone se situait globalement entre les 


cotes 490 et 450 NGM (voir annexe 1.f). 


2: 


Géologie et géotechnique : 


Le laboratoire public d’essais et d’études (LPEE) s’est mis sur place pour effectuer des essais 


in situ afin d’aboutir aux logs stratigraphiques, les profils lithologiques (voir annexe 1.g et 


1.h), et les procès-verbaux des essais in-situ. Sans oublier de compléter l’étude par des essais 


en laboratoire. En effet, les matériaux meubles et rocheux prélevés à partir des cavaliers ont 


été soumis aux essais de laboratoire suivants : 


V 


SSKK SOS SSS 


V4 


Granulométrie par tamisage. 

Essai de sédimentométrie 

Double sédimetométrie 

Limites d’atterberg 

Valeur au bleu de méthylène 

Essai mécanique sur sol (essai triaxial) 

Essais mécaniques sur roches de type Los Angeles et Micro Deval 
Essai Proctor 


Essai Pin hole test 


$Æ Matière organiques dans les cavaliers : 


Compte tenu des teneurs quasi-nulles en matières organiques (0% < MO < 3 %), on conclut 


que les remblais analysés sont inorganiques. 


$Æ Le site de la future digue et son périmètre environnant montre une relative 


homogénéité dans la répartition latérale et verticale des faciès traversés. 
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Æ La colonne lithologique est marquée par la présence de phosphates gréseux 
entrecoupés de bancs de calcaires et par une couche de 0.6 à 2.50 m de marne 
phosphaté jaunâtre. 

Æ Les matériaux des cavaliers sont abondants et majoritairement classées C1B5, C1B6. 
Ces sols sont généralement réutilisables en remblai selon la classification GTR (Guide 


des Tracés Routiers). 


Sables limono-graveleux 
avec des gros éléments 
Limons sablo- graveleux 






avec des gros éléments 


Tableau I: types de matériaux des cavaliers 


> Donc on aura une quantité suffisante des matériaux de construction de bonne 


qualité à partir des cavaliers. 


HI. Choix du type d'ouvrage : 


Les variantes envisageables pour la future digue : 

v Le remblai homogène : dont le corps principal est constitué d’un seul matériau avec 
des protections en petit enrochement en aval. L’étanchéité reposera soit sur le corps de 
la digue (matériau argileux ou limoneux), ou bien sur d’une couche de géo-membrane 
posée sur géotextile pour la protéger contre le poinçonnement. Le complexe sera ancré 
en crête et au pied du talus amont de la digue. 

v Le remblai zoné : comportant un noyau étanche (épais ou mince) en matériaux 
argileux ou limoneux, enveloppé par des transitions (fltre/drain), avec des recharges 
amont et aval en alluvions tout venant. 

En se basant sur l’estimation des quantités de fines existantes sur terrain et sur l’étude 
comparative de cout entre les deux variantes, 1l s'est avéré raisonnable de retenir un profil de 
la digue comportant : 

v Un remblai compacté sans tri de matériaux, muni d’un complexe d’étanchéité 
(géomembrane + géotextile) (voir annexe 1.c et 1.d). En plus de sa fonction 
protectrice, le géotextile jouera le rôle de filtre synthétique. Un drain-cheminée en 


matériau drainant extrait de la zone des cavaliers par criblage (élimination des fines). 
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Cette coupe type est valable sur la partie sud de la digue, où sera implanté le 


système de récupération des eaux. 


AXE DE LA DIGUE 










GECMEMBRANE+GEGTEXTILE 


REMBLAIS COMPACTE 


280 CLARINETTE 
CONOUITE METALLIQUE 
7 ES mm 


| EXTERIEUR | 2.80 


DRAIN 


CONDUITE EN ACIER NOIR & 500m POSEE 
SUR DES SOCLES EN BETON 


N a 
R- y 


ÉTRTITEER EIRE OIPPIEET EANET EISELI TION TELTI ENEE EIT EIRUE NEEESE PPELAIEELE EDIEN ETEENI UEN ; 
na BETON D'ENROBAGE 


Figure 6: coupe type de la digue au Sud du bassin 


La géomembrane offre l’avantage d'empêcher l’infiltration des eaux à travers les remblais. Il 
est nécessaire de prévoir une évacuation immédiate de toute eau filtrée afin de réduire la 
percolation à travers les remblais et par conséquent, éviter le développement des pressions 
interstitielles dedans (risque d’instabilité). 
v Un remblai zoné comportant les matériaux à taux élevé de fines (>35%) sur les 
premiers cinq mètres. Les remblais à taux de fines réduits seront mis en place sur le 


reste du corps de la digue. 


Cette disposition sera valable sur l’ensemble de la partie centrale et nord de la digue. Les 
remblais étanches remplaceront la géo-membrane et seront en contact avec les boues en 


dépôt. 


| INTERIEUR | 


13 m 





Figure 7: coupe type de la partie centrale et nord de la digue 





LE, A 


Mémoire du projet de fin d’études —/ AN 


ECOLE MOHAMMADIA D'INGENIEURS 


IV. Mécanisme d’alimentation de la boue et récupération des eaux : 


La forme du bassin sera rectangulaire pour prolonger le chemin d'écoulement de la 
boue et des eaux, favorisant ainsi la décantation des fines et la récupération d’une eau 
suffisamment claire à l’autre bout du bassin. 

En s’inspirant de la technique adaptée sur les digues d’El Merah (MEA), notre bassin sera 
alimenté par 3 points implantés sur le mur en remblais, la distance entre ces 3 points va être 
abordée dans la partie simulation sur RIFT TD afin de trouver la bonne disposition qui nous 
donne une répartition uniforme de la boue. 
Quant à la récupération, trois systèmes sont également identifiés sur les digues, à savoir : 

v Cheminées drainantes (constituées d'un ensemble de clarinettes) (voir annexe 1.e) 

v Pompes flottantes (au cas où il y'a de l'eau stagnante) 

v Remblais drainant/ filtrant 
Pour la présente digue de Benguerir, le choix est porté sur une combinaison des deux premiers 
systèmes, avec une amélioration des cheminées consistant à les mettre sur le parement amont 
de la digue, moyennant une assise en béton (pour l’aménagement de marches d’escalier). Il 


y’aura 4 points de récupération des eaux. 








GARDE-COAPS 
METTALIQUE 


ESCALIERS D'ACCES 
AUX CLARINETTES 


Figure 8: escaliers d'accès aux clarinettes 


Les clarinettes seront connectées à leurs bases à une conduite collectrice logée dans une 
tranchée et enrobée dans le béton, afin d’éviter sa détérioration en raison du tassement des 
remblais. (Voir figure 12). La conduite collectrice a une pente de 0,9% afin de favoriser 


l’écoulement de l’eau vers l’extérieur du bassin. 
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V. Traitement de la fondation : 
D’après les essais (Lugeon et Lefranc) réalisés par LPEE sur la fondation, il s’avère 
que le terrain en total est imperméable (valeur du coefficient de perméabilité est inférieur < 
SUL), avec nécessité de traitement par injection non couteux dans quelques zones vu les 


valeurs proches aux limites de l’imperméabilité. (Voir annexe 1.1) 


VI. Choix du dispositif d’auscultation : 


En plus de l’inspection visuelle systématique de la digue et son environnement (la 
crête, les talus et le pied de remblai), consistant en un repérage fuites, suintement, tassement 
ou fissures, La conception de la digue doit comporter un dispositif d’auscultation. 

Le but d’un dispositif d’auscultation est de mesurer des phénomènes ayant une incidence 
majeure sur le comportement de la digue. 
Par sa définition, l'auscultation est constituée d'un ensemble de données quantifiées, fournies 
par des instruments de mesure, qui permettent de suivre le comportement de l'ouvrage. 
Le système d'auscultation de la digue doit être conçu de manière à surveiller les paramètres 
suivants : 
e Les déformations des remblais de la digue pendant la construction et à 
l'exploitation. 
e Les pressions hydrostatiques dans la fondation de la digue. 
e Les pressions interstitielles qui se développent au contact du noyau 
(éventuellement) avec la fondation. 
e Les fuites collectées à l'aval de la digue. 
Pour ce faire, on prévoit l’installation des piézomètres le long de la digue afin de maitriser les 


pressions interstitielles développées à cause du contact direct de l’eau avec les cavaliers. 
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Figure 9: la disposition des piézomètres le long de la digue 
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Figure 10: coupe a-a type d'un piézometre 
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Chapitre III : Caractérisation de la boue de lavage de 
phosphate à la sortie de la laverie 
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A la sortie de la laverie, la boue est dirigée vers un épaississeur par gravité pour 


récupérer une quantité importante de l’eau claire. 


Le laboratoire américain RADISE s’est chargé des essais d’identification des caractéristiques 
de la boue à l’entrée et à la sortie du décanteur. On va s’intéresser à la composition de la boue 


à la sortie du décanteur (après épaississement). 


I. Analyse granulométrique : 


L’analyse granulométrique permet de déterminer et d’observer les différents diamètres de 
grains qui constituent un granulat. Pour cela l’analyse consiste à séparer et classer à l’aide de 
tamis ces grains selon leur diamètre. Les grains ainsi isolés de l’analyse permet d’observer et 
d’exploiter ces informations très simplement. Les manipulations et les conditions de 
manipulation sont décrites par la norme NF P 18-560. Elle concerne les granulats d’un 


diamètre supérieur à 100um. 
Ainsi le tableau 1 représente les résultats de l’analyse granulométrique par tamisage. 


SIEVE ANALYSIS 


No. 4 
No. 10 
No. 20 
No. 40 
No. 50 
No. 60 
No. 100 


No. 140 


No. 200 


Hydrometer Analysis 





Tableau 2: analyse granulométrique par tamisage 


À partir des essais de tamisage tableau 2, on trace la courbe qui représente le cumul passant 


en fonction de diamètre des grains sur une échelle logarithmique, on constate à partir de la 
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courbe que les particules de diamètre inférieur ou égal à 250um constituent un pourcentage de 


95,6%, c’est-à-dire après lavage des phosphates, on coupe à peu près à 200um. 


cumul passant % 


IL. 


100 
90 
80 
70 
60 
50 
40 
30 
20 
10 


courbe granulométrique après épaississement 








250um; 95,6% 


74um; 54% 


1 0,1 0,01 0,001 


Figure 11: courbe granulométrique de la boue 


Vitesse de chute : 


La vitesse de chute est calculée à partir de la théorie de Kynch : 


Ac 


Qu : 


Cu 
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Area for — Co | C, 
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Figure 12: Paramètres caractérisant le décanteur 


: débit entrant au décanteur (m3/s) 
Co : 
Ce : 


concentration de solide suspendue à l’entrée kg/m? 


concentration de solide suspendue à la sortie en haut kg/m’ 


- surface de clarification en m? 


: surface de décantation en m? 


débit de boue sortant du décanteur (m°/s) 


: concentration de solide décantée sortante kg/m? 





ES) A 
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Vi: vitesse de chute après épaississement (m/s) 
Ci: concentration de solide dans la couche de décantation kg/m° 


Ub: vitesse de pompage de la boue à la sortie (m/s) 


X C 
= . /c x A, — Up 


En utilisant cette relation, on développe le tableau 3 suivant qui nous donne la vitesse de 


chute dans l’épaississeur en fonction du diamètre des grains. 


SETTLING VELOCITIES 


1653.0 
291.8 
31.5 
12.9 
6.4 
4,5 
1.62 
0.80 
0.40 
0.0 747 
0.0322 
0.0047 
0.0026 
0.0012 
0.0007 
0.00007 


< 0.00007 





Tableau 3: vitesse de chute en fonction du diamètre des grains 


A partir du tableau 2, on trace la courbe représentante la vitesse de chute en fonction du 
diamètre des grains sur une échelle logarithmique, on constate que la vitesse de chute peut 


arriver à 2 m/s pour la portion 0,1% de diamètre égal à 2mm. 
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vitesse de chute après épaississement en (cm/s) 
10000 


1000 









o 4,76 mm; 1653 cm/s 
10 
1 
0,1 0,074mm; 0,4cm/s 
0,01 
0,001 
0,0001 


0,00001 
0 0,5 1 1,5 2 2,5 3 3,5 4 4,5 5 


dimètre des grains en (mm) 


Figure 13: courbe de la vitesse de chute en fonction du diamètre des grains 


IHI. Caractéristiques du produit boue : 

1. Teneur en eau 
On désigne par teneur en eau (moisture content en anglais) la quantité d'eau liquide contenue 
dans un échantillon de matière, cette quantité étant évaluée par un rapport pondéral au sens 


ingénieur ou volumétrique au sens physicien. Elle est exprimée par : 


Wọ, : Poids de l’eau en Kg 


W; : Poids des grains en Kg 


w : Teneur en eau exprimée en % 


Moisture Content (%) 





2. Gravité spécifique 
La gravité spécifique est le rapport entre la masse volumique d'une substance (partie solide) et 
la masse volumique d'une substance de référence (eau). De manière équivalente, 1l s'agit du 
rapport entre la masse d'une substance et la masse d'une substance de référence pour un même 


volume donné. 
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ern Psolide 





PH20 


Specific Gravity 





3. Limites d’Atterberg : 


Le but de l’essai des limites d’Atterberg est d'identifier et classifier les l’échantillon, et pour 
déterminer La limite de plasticité on utilise la "Méthode du rouleau». Quant à la limite de 


liquidité, se détermine par la "Méthode de la Coupelle de Casagrande". 





1. Appareillage : 


y 





Figure 14: appareillage de détermination des limites d’Atterberg 


2. Mode opératoire : 


Cet essai est réalisé sur la partie granulométrique appelée mortier inférieur à 0,4mm 
e On malaxe rigoureusement la totalité de l’échantillon afin de bien homogénéiser 
e On remplit la coupelle au tiers. 
e On trace le milieu de l’échantillon à l’aide d’un outil à rainurer et l’essai comme 
immédiatement 
e La coupelle fixée à l’appareil est soumise à une série de chocs réguliers jusqu’à ce que 
les deux lèvres se ferment (le nombre de coups de fermeture doit être compris entre 15 
et 35) 
Pour confirmer le nombre de chocs, il faut recommencer immédiatement l’essai, si les deux 
essais successif ne différent pas plus d’un choc, on prélève à l’aide d’une spatule deux 


échantillon de chaque côté des lèvres et on détermine la teneur en eau. 


3. Calcul des limites d’Atterberg : 


W. : la limite de liquidité est la teneur en eau correspondant à la fermeture à 25 coups. 
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Pour déterminer Wp, on prend une partie du mélange et on essaye de faire des petits 
rouleaux sur une plaquette jusqu’à obtenir un diamètre de 3mm. 

Après on soulève le rouleau à 1 ou 2cm au-dessus du plancher, quand il se fissure et 
rompt. On prélève sur chaque morceau des échantillons pour déterminer la teneur en eau qui 


sera la valeur de Wp, Et il en vient l’indice de plasticité Ip = Wi - Wp 





Pour le cas de notre produit boue, on trouve : 


Atterburg Limits 
Liquid Limit (%) 


Plastic Limit (%) 


Plasticity Index 





4, Classification USCS : 
La classification USCS (Unified Soil Classification System) (1965) utilise les résultats 


d’essais classiques d’identification des sols : 


e des critères granulométriques : les pourcentages de gravier, sable et particules fines 
(tamisats à 2 mm et 0,08 mm); la forme de la courbe granulométrique : coefficient 
d’uniformité ou de Hazen Cu, coefficient de courbure Cc (ou Cz). 

e les caractéristiques de plasticité Wx et Ip, et la ligne A d’équation : 


Ip = 0,73*(Wz — 20) (relation de Casagrande). 


e la teneur en matières organiques. 


La classification peut également s’effectuer à partir de l’observation visuelle du sol et de tests 
simples de chantier. Mais il faut une grande expérience pour appliquer correctement cette 


méthode de classification de chantier. 


Æ Les groupes de sols : Le système de classification USCS débouche sur 15 sols types, 


affectés chacun d’un symbole à deux lettres, prises dans les trois ensembles suivants : 
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p Granularité du | 
Eléments du sol i Plasticitė du sol 


G : Grave. Le gravier est la fraction 


principale 

S : Sable. Le sable est la fraction 

principale b : Bien gradué t : Très plastique 
L : Limon ou limoneux m : Mal gradué p : Peu plastique 
A : Argile ou argileux 

T : Tourbe 

O : Organique. Le sol contient des 


matières organiques 


Tableau 4: groupes de sols 


+ La procédure de classification : La figure 15 donne la classification des sols fins et 


le tableau 5 la classification des sols grenus. 








Figure 15: classification des sols fins 
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Définitions 


moins de 
5% 
Plus de 
50% des 
< 0,08 mm 
éléments > 


0,08 mm 


d'éléments 


ont un 


| plus de 
diamètre > 


12% 


< 0,08 mm 


amm 


Plus de 
50% des 5% 


éléments > | d'éléments 


moins de 


plus de 
1290 
d'éléments 


< 0,08 mm 


Cu = Déo/ Dio > 4 et 
1< Ce = (D:0}°/D10D60 
< 3 
Une des conditions 
de Gb n'est pas 
satisfaite 
Limites d’Atterberg 
au-dessous de la 
ligne A 
Limites d’Atterberg 
au-dessus de la ligne 
A 
Cu = Déo/Di0 > 6 et Sable 
1< Ce = (D:0)°/D:0D4 | propre bien 
< 3 gradué 


Une des conditions 
de Sb n'est pas 
satisfaite 
Limites d'Atterberg 
au-dessous de la niti 
ligne A 
Limites d’Atterberg 
au-dessus de la ligne 
A 


limoneux 


Lorsque la teneur en particules fines (< 0,08 mm) est comprise entre 5% et 
12% on utilise un double symbole. Par exemple : Sb-SL 





Tableau 5: classification des sols grenus 


Pour notre cas, on a W:=46,7% et Ip=23,1% donc selon la figure 15 la boue est considérée 


comme argile peu plastique et on a selon la figure 11 : plus de 50% des éléments >0,08mm et 


ont un diamètre < 2mm, donc selon le tableau 5 la boue est un sable et puisque plus de 12% 
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d’élèments < 0,08 mm selon la figure 11, donc on a le choix entre les classifications : SM et 


SC. 


Et on selon la figure 15 le point représentant de la boue est au-dessus de la ligne A, d’où la 


classification de la boue selon la norme USCS : 


USCS Classification 





Clayey sand (SC) 


5. Composition du total solide : 


Le total solide est déterminé à partir de la teneur en eau par : 


TS = : 
= 1+w 

Puisque : 

Wy Wr — W; 1 

Ww = — = = Å= — 1 

Ws Ws TS 
Avec TS : total solide en % 
w: Teneur en eau en % W; : Poids solide en Kg 
Wy : Poids de l’eau en Kg Wr : Poids total de l’échantillon (sohide+eau) en Kg 


Pour notre produit boue, on trouve : 


Total Solids (%) 





Remarque : on constate que les essais faits par le laboratoire externe RADISE nous donnent 
un pourcentage de solide de 37,38% après épaississement, or le bureau d’études JACOBS S.A 
a dimensionné le décanteur et la digue d’épandage à un pourcentage de solide de 35%, qu’on 


va utiliser dans la simulation. 
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Chapitre IV : Modélisation du comportement de la boue 
à l’intérieur de la digue 
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Si le modèle de fluide newtonien décrit bien la très grande majorité des fluides 
composés de molécules simples, 1l existe un bon nombre de fluides, dont certains sont d'usage 
très courant, qui ont un comportement sous l’écoulement plus complexe. La définition d’un 
fluide newtonien est assez restrictive : les contraintes de cisaillement sont proportionnelles au 


gradient de vitesse, ce qui implique que : 


e Dans un écoulement de cisaillement simple, les seules contraintes crées par 
l’écoulement sont des contraintes de cisaillement. 

e La viscosité est indépendante de la vitesse de cisaillement. 

e La viscosité est indépendante du temps et les contraintes s’annulent immédiatement 


lorsque l’écoulement est arrête. 


Toute déviation de ces règles est le signe d’un comportement non-newtonien. La 


description de ces comportements et leur interprétation en relation avec la structure 


microscopique du fluide constitue la discipline appelée Rhéologie. 





I. Hypothèses de l’étude : 


œGravel| %Sand %silt Clay 





(i 46 21.5 32,5 
Tableau 6: composition de la boue 


e [a boue avec un taux de solide de 35% est considérée comme fluide non newtonien 
rhéo-épaississant (shear thickening). 
e Avec un taux d’argile (Clay) égal à 32.5% (> 10%), la boue est considérée comme 


matériau cohésif. 
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II. Etude théorique de l’écoulement de la boue à l’intérieur de la digue : 


Pour ce faire, on va diviser le bassin en trois parties : 


1°" étape : Jet de la boue 
dans le bassin 


2°"e étape : 


3°M€ étape : écoulement de 
l'eau vers les points de 
collecte. 


e érosion et dépôt du 
solide. 
Décantation, 
consolidation et 
drainage. 





1. Etude de l’écoulement de la boue à l’entrée de la digue (théorie du JET) : 


Pour le cas de la digue d’épandage de Benguerir, les conduites d’entrée utilisées sont 


circulaires de diamètre 400 mm, donc on va étudier le cas d’un jet circulaire. 


Le jet est bidimensionnel puisque les conduites vont être allongées le long du talus nord du 


bassin. 
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ÉYIRTUAL 
ORIGIN : 


Figure 16: zones de développement d'un Jet 


e ‘flow development region ‘ connue également par la zone de développement du jet, 


sa longueur est estimée de [0 ; 6-10Do]. Pour notre cas ça donne [0 ; 4m]. 


e ‘fully developped flow region’ connue par la région du jet totalement développée, 


elle s’étend de 4m à l’infini. 


La vitesse et la concentration du taux de solide suivant l’axe du Jet : 








U . D 
“DE —_6 9-0 
Ug X TE 
in D = + 
C nax D, Pki 
Co X 


La vitesse des particules à une distance r de l’axe du Jet et à une distance x de l’entrée : 


= exp Ene 


u 








QUE 


| a 

c ro] 

-E exp|—62 s 
X 


Cmax 


v La vitesse d’entrainement radial v ne dépasse pas les 5% de la vitesse axiale maximale 


Umax 


=i lergih unii 
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Ces équations sont développées à partir des équations de Reynolds exprimées dans le système 


cylindrique. 


2. Applications des équations du jet pour notre cas : 


Pour Do=400mm, u0= 1,65 m/s, Co= 1,31 Kg/L, on a : 


vitesse axiale du jet 


2,5m ; 1,61 Kg/L 


O 
`~ 
5 
x 
© 
€ 
=) 


13,4 m; 0,3 m/s 


10 15 


distance de l'entrée en m 





Figure 17: vitesse de la boue suivant l'axe du jet 


On constate que la vitesse maximale suivant l'axe reste à une valeur de uo sur une distance de 


2,5m, après elle diminue pour arriver à 0,3m/s à une distance de 13,4m. 


concentration axiale du solide 


2,5 m; 1,13 Kg/L 


Cmax en Kg/L 


10 15 


distance x en m 





Figure 18: concentration du solide suivant l'axe du jet 


On constate que la concentration en solide diminue rapidement à une valeur de 142g/L à une 


distance de 20m à partir de l'entrée. 
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1,8 
1,6 
1,4 
12 


vitesse axiale aux 
1 bords du jet à x=2,5m 
0,8 


0,6 concentration du 
0,4 solide aux bords du 
jet à x=2,5m 
Le r= 0,6m 
O 
-0,2 0 0,5 il il 
distance r prp à l'axe du jet en m 





Figure 19: caractéristiques du jet à 2,5m de l'entrée 


On constate que la vitesse et la concentration en solide diminuent rapidement pour s'annuler à 
une distance de 60cm à partir de l'axe du jet. Donc le jet est encore concentré autour de son 


axe principal. 


vitesse axiale aux 
bords du jet à x=6m 


concentration du 
solide aux bords du 
jet à x=6m 





Figure 20: caractéristiques du jet à 6m de l'entrée 


On constate qu'à une distance de 6m à partir de l'entrée le jet commence à s'élargir et avoir 


une distribution de solide plus étalée. 





KA 
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1. Développement des équations de Saint Venant: 
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Dans le cas bidimensionnel, pour décrire la hauteur de l’eau dans le bassin on utilise les 


équations de Saint-Venant : 














öh , À r | à . | E 
Se tpr (E +h) + z (+ hu) =0, 

Ou ðu + „2! f öh. bu 4 (a E Lu 
g TU TA = UT 

ot * Oz y Tor Ga D 
dv Ov gu | Fy Fy 
& "os ‘5 + fu 58 -wr (Et + 2 


AVES : 


ÉD DT <= S 


Vitesse suivant x, direction de l’écoulement 

Vitesse suivant y 

Hauteur entre la surface libre d’eau et le terrain 

EÉlévation de la surface du terrain par rapport à la référence 
L’accélération due à la gravité 

Coefficient de Coriolis associé à la force de Coriolis 
coefficient de la trainée visqueuse 


La viscosité cinématique 





Figure 21: schéma des paramètres utilisés dans les équations 
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2. Développement de l’équation d’Exner pour la dispersion des oui. : 
L’équation qui donne la variation du lit n en fonction du temps t, du coefficient de dispersion 


£o et du débit solide qs. 


IN 1 
=- VẸ.re 
dt Ep E 


1 _ ppe | 
AVEC : Ep = E tel que pọ est la porosité des sédiments du lit. 
7 FO 


L'étude du comportement de la boue à l’intérieur du bassin est obtenue par le 


couplement des équations de Saint-Venant et équation d’Exner (voir la figure 21 pour 





la notation utilisée): 























ðh. 
ot aa =e 
Ju Maia i  — — fm — 
I "Dz Dy - 9% 
Où nu. du om S FN 
& “Drz Əy TD ðr? Əy? } ` 
dn l — —> 
E — => — * Cle 
Ot Eao 


3. Caractérisation de l’écoulement à l’intérieur de la digue : 





On a trouvé que la vitesse après 6 m du nord du bassin (zone du jet totalement développé) ne 


dépasse pas 0,7m/s, on a le nombre de Froude suivant : 





Fr =- 7016 <1 m Régime fluvial 


Avec : U=0,7m/s, D=4*Rp, g=9,81m/s?, = 1076m? /s , h=0,5m (on a intérêt à garder la 


hauteur d'eau faible). 


= | b+m-h)-h 
NA n AL 


« - — b+2.-h\1+ m? 











Puisque b est très supérieur à h, on estime le rayon hydraulique Ra à h. 


On calcule le nombre de Reynolds : 





Mémoire du projet de fin d’études 


LL U*D 


R = ~ 1400000 > 2000 = écoulement turbulent 





4. Transport de sédiments en régime fluvial : 
1. Classification du transport de sédiments : 


On distingue deux types : 


Wash load (transport par lessivage) : correspond aux particules très fines (silt et argiles) qui 


demeurent en suspension quel que soit le régime d'écoulement fluide. 


Bed material load (transport du fond et des berges) : correspond au transport de sédiments 


(sable et gravier) qui demeure lié aux matériaux du lit et aux berges et se manifeste en 


transport durant les écoulements turbulents. 


=} Modes de transport dans un écoulement turbulent : 





Charriage de fond (roulement et glissement) 
Saltation (sauts successifs) 


Suspension (transport permanent au sein du fluide) 


KO OK OS 


Solution (transport chimique) 





ue, Rolling 











Bedload 
td Saltation 
A aenn ARE ee san, Te ig, Suspension Suspended 
load 





Figure 22: différents modes de transport des sédiments 


2. Interprétation hydraulique de la distribution granulométrique : 
Diagramme de Shields : 
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Le mouvement d’une particule se produira quand les forces instantanées exercées par le fluide 
sur la particule sont juste supérieures aux forces de résistance instantanée apparentée au poids 
submergé de la particule et au coefficient de frottement. Les forces d’entraînement sont 
fortement liées aux vitesses locales près du fond. Dans des conditions turbulentes de 
l’écoulement, les vitesses sont fluctuantes dans l’espace et dans le temps, ce qui montre, avec 
l’aspect aléatoire de la taille des particules, la forme et la position, que le début de mouvement 
et d'entraînement n’est pas seulement un phénomène déterministe mais aussi un processus 


stochastique. 


Bonnefille(1963) et Yalin (1972) ont montré que la courbe de Shields peut être exprimée 


d’une façon plus simple, en introduisant le diamètre adimensionnel de la particule D, : 


(s-1)* g"? 
Avec : D:: diamètre adimensionnel 


s : gravité spécifique adimensionnel 

V : viscosité cinématique de l’eau en m?/s 
D5o : diamètre à 50% de passants en m 
g=9,81 m/s? 


Selon la courbe granulométrique figure 17 on a : 
Dso = 0,074 mm ; s= 1,31 ; 9=10% m°/s; g=9,81 m/s? on trouve: D:=75 (adimensionnel) 


Ainsi, la courbe de Shields peut être représentée par : 


44 
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Figure 23: Début de mouvement et de suspension sur un fond plat 
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u; X m . 
Avec : Oer = a paramètre critique de Shields 
La contrainte de cisaillement au fond : Ther = p *UËr 


En se basant sur D:calculé et à partir du diagramme de Shields (courbe de Van Rijn), on 


trouve la valeur critique de début de mouvement: 


Log (D:)= 1,87 ; la courbe nous donne donc : log(8)=1 donc Ocr=2,718 
Ucr = ,/@cr * (s — 1) x g * dso = 0,025 m/s 


cette valeur reflète une idée sur la granulometrie fine de la boue, et puisque la vitesse à 
l’entrée est de l’ordre de 0,7 m/s et l’écoulement est turbulent donc on aura une érosion forte à 


l entrée. 


3. Dépôt: 
Le dépôt est prédominant quand la contrainte de frottement au fond ( tb ) est inférieure à la 
valeur critique de dépôt( td ), comme 1l a été montré par Krone (1962 ). Il a étudié le processus 


de dépôt des matériaux cohésifs en appliquant différentes vitesses générant une circulation de 
la vase naturelle dans un canal du laboratoire. 60% des particules composant le matériau (vase 


de San Francisco Bay) avaient une taille inférieure à 2 um, et 50% avaient une taille inférieure 





EL? 
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à | um. Krone a montré que le taux de deposition du solide est fortement proportionnnel à la 
concentration en sediments. 


107? 
| 





settling velocity (m/s) 
deposition rale (kg/ms) 


o 


hindered 
settling 


1072 107! 109 0! 10? 
——+ sediment concentration (kg/m?) 


Figure 24: Taux de dépôt en fonction de la concentration 


4. Erosion : 


Les particules sédimentaires, les floques ou les morceaux de la surface du fond subiront des 





érosions quand la contrainte de cisaillement au fond appliquée (tt ) dépasse la valeur critique 


d’érosion( Te ), qui dépend des caractéristiques du matériau du fond( composition minérale, 
matériaux organiques, salinité, densité etc. ) et de la structure du fond. Des résultats 
expérimentaux et des recherches ont montré que cette valeur critique est supérieur à la valeur 
critique de dépôt( Te > ta ). Mehta et Partheniades( 1975 ) ont constaté que même pour des 
contraintes de cisaillement supérieures à 1.5 N/m? dans un écoulement stationnaire, aucune 
érosion n’est observée, car le fond a été composé de floques déposés de forte densité. Ainsi, te 
peut être très grand dans le cas des couches supérieures du fond composées de forts dépôts de 


floques. 


Plusieurs types d’érosion sont distingués dans la littérature : 
1. L’érosion de surface qui correspond à une érosion de particules une à une ou sous forme 
de floques. 
2. L’érosion de masse qui consiste en mouvement (pick-up) de floques ou de fragments à 
cause de fissures dans le fond. Plusieurs tentatives ont essayé de lier te à certains 
paramètres de base comme l’indice de plasticité, l’indice des vides, la teneur en eau, la 


cohésion etc. 


KA 
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1. Erosion des matériaux récemment consolidés : 
Pour les matériaux récemment consolidés, la contrainte d'érosion critique est calculée à partir 


de la vitesse critique, on trouve: 


Z 2 
Te = P * Ucr 


p =10° Kg/n? ; u,, = 0,025m/s 


0,63N /m? 


Te 


2. Erosion des matériaux très consolidés : 
Pour les matériaux très consolidés, on parle de probabilité d'érosion qui dépend de la vitesse 
de l'eau au point considéré et la hauteur d'eau au-dessus du point, par exemple, pour une 
vitesse de 0,7 m/s (trouvée à une distance de 6m de l'entrée) et une hauteur d'eau de 8cm, on 


est à 50% sûr que les particules très consolidées vont être érodées. 


CLAY SOILS 


Lu d 
Les 


= 
> 


probability of channel 





Cas J a) nm Un 


CHANNEL VELOCITY ft/sec 


-- stable 


o-seoure 


À 2 3 45 l 2 345 10 20 
DEPTN OFFLOW ft 


Figure 25: Vitesses critiques d’érosion pour des dépôts argileux très consolidés, Hughes(1980) 


5. Écoulement de l'eau vers les points de collecte: 
En utilisant la formule de Manning-Strickler donnant la vitesse moyenne de l'écoulement, on 
calcule la vitesse pour différentes pentes, afin de pouvoir comparer le temps nécessaire pour 
qu'une goutte d'eau arrive à la sortie après drainage de la boue : 


x į¿1/2 


V=K*R 


> WIN 
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Nature des parois Valeur de K en m?” /s 


CS 


Rivière à berges étroites très végétalisées 
Lit majeur en prairie 

Lit majeur en vigne ou taillis 

Lit majeur urbanisé 


Lit majeur en forêt 





Tableau 7 : coefficient de Manning-Strickler 


On prend : 
K= 13 m/s 





i = 1%o V = 0,26m/s 

Plus la pente est grande plus la 
i = 1% V = 0,82m/s récupération d'eau est plus 
i = 2% V = 1,16m/s rapide. 


MI. Simulation sur RIFT TD 

1. Présentation du logiciel RIFT TD : 
Rift TD est un modèle numérique de terrain tridimensionnel avancé spécialement conçu pour 
modéliser le dépôt de résidus. C'est un outil puissant qui facilite l'installation de stockage des 


résidus, 1l permet : 


e évaluation des options / sélection du site 
e conception conceptuelle et détaillée 

e planification des dépôts 

e fonctionnement et gestion 


e Gestion de l'eau 


Rift TD est utilisé souvent pour modéliser la déposition des résidus à l’amont et à l’aval des 
barrages et la gestion des déchets dans les bassins de rétention tout en modélisant les 
variations des volumes d’eaux, c’est un logiciel payant mais la version complète existe 


également sous licence d’évaluation. 
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V Fonctionnement du logiciel : 


Après introduction des données nécessaire pour lancer la simulation, on choisit un pas 


d’élévation entre la déposition n et la déposition n+1. Rift TD lance le jet des résidus et au 


cours de la déposition n par exemple, les courbes de niveaux caractérisent l’élévation de la 


boue jusqu’à atteindre l’élévation de la déposition n+1 et le logiciel nous donne les 


informations propres à la déposition n (la durée de la déposition, le volume et le tonnage 


cumulé, le volume de l’eau généré..). 


2 


Présentation des données : 

Le débit solide : 4573 tonnes/jour (190,51 tonnes/heure) 
Le débit de l’eau : 544,31 tonnes/heure 

Le débit de boue : 734,82 tonnes/heure 


Ce sont les débits nominaux de dimensionnement des conduites liant le décanteur à la digue 


d'épandage. 


La densité après consolidation et drainage (55% de solide) du profil subaérien : 





1,4 t/m? (sub-aerial dry density) 


Selon le bureau d’études Team Maroc, le but est d’atteindre une consolidation à 55% de 


solide sur les trois ans premiers pour la digue de Benguerir, qui correspond à une gravité 


spécifique de 1,4 t/m°. 


La densité après consolidation et drainage du profil subaquatique: 1,2 t/m? (sub- 


aqueous dry density) 


Profil subaëérien 


Profil subaquatique 


ER 


Eau récupérée 





Figure 26: profil subaérien et subaquatique 
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e La pente choisie pour la plage subaérienne est fixée par expérience à 5%o pour ne pas 
avoir un tassement différentiel lors de la consolidation, et celle de la plage 
subaquatique est égale à 2% afin d’avoir un drainage rapide de l’eau. 

e Le volume d’eau récupéré est maitrisé par la revanche (freeboard) de stabilité estimée 


à 2m. 


Les points d’entrées varient suivant la pente du talus (1 V/2,6H) verticalement de 472NGM à 


485NGM (la hauteur du bassin= 13m =485NGM - 472NGM), le pas d’élévation dépend de 


l’incrémentation choisie, on prend 50 incrémentations c’est à dire : 


Pas =13/50=26cm équivalent à 73cm suivant le talus. 





Tableau S: données utilisées pour la simulation 


Le début de simulation choisi : 18/04/2019 


Les scenarios étudiés : 
nicane 
| 
sans chicanes 


sans chicanes 
avec avec chicanes — 


sans anaa 


avec chicanes 
sans chicanes 
avec chicanes 


Tableau 9: scénarios étudiés 





3. Présentation des résultats : 


1. Scenario 1 : Digue d’épandage à pente 0% sans chicanes 


vV Géométrie : 


Mémoire du projet de fin d'études =/ AN 
LE, 


A 


ECOLE MOHAMMADIA D'INGÉNIEURS 





Figure 27: géométrie de la digue à une pente 0% 


En utilisant VBA excel, j'ai pu programmer un algorithme afin de dessin le bassin pour toute 


pente voulue avec l’option d’ajouter un nombre de chicanes. 





0 Sans 68 41,33 3 ans 7 mois 20 jours 


Tableau 10: résultats de récupération d'eau de la digue à une pente 0% 


vV Les formules de calcul : 
La durée de la 1%° récupération est donnée par le logiciel. 


D volume collecté en n jours (m°) 
o TÉCUPÊTÉ par n jours = ———— LLa 
pere p J n + débit de boue journalier 


La durée de saturation est donnée par le logiciel. 


On constate que cette configuration du bassin n’est pas bénéfique vu que le taux de 


récupération ne dépasse pas 42%. Concernant la durée de la 1% récupération, on constate 
qu’elle dépasse 2mois, c’est normal vu la longueur du bassin qui dépasse 1km et l’absence 


d’une pente accélérant l’écoulement vers les points de collecte. 


V Le bassin après saturation : 





Mémoire du projet de fin d'études =/ W 


pus) A 


ECOLE MOHAMMADIA D'INGÉNIEURS 
DTM View  3D View Results Verification Results 
ny X 
Date: 08 déc. 2022 
Raise Volume: 144688 
Raise Tonnage: 186 160 
Cumulative Volume: 4625 347 
Cumulative Tonnage: 6 084 793 





nf ESRI; 


Figure 28: l'élévation de la boue après saturation de la digue à une pente 0% 
Date: 08 déc. 2022 
Raise Volume: 144 688 
Raise Tonnage: 186 160 


Cumulative Volume: 4 625 347 








| 
Tableau 11: taux de stockage de la digue à une pente 0% 


volume de solide stocké m° 


% de stockage = ——————————————————————— 
j g volume total du bassin m? 


On constate qu'avec un pourcentage de 52% de stockage de solide (tailings) on perd presque 


la moitié en terme de capacité de stockage du bassin. Donc ce scenario est à éviter. 


V L'élévation du bassin en longueur après saturation : 


Mémoire du projet de fin d’études =/ AN 


LEE}, A 
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Figure 29: L'élévation du bassin en longueur après saturation 


Remarque : le chainage veut dire la distance en (m) 


On constate que la boue consolidée et drainée n’est pas bien répartie, en fait, la boue se 
décante sur les 600 m premiers, d’où la grande perte en terme de stockage. Un autre problème 
concernant la pente raide de la boue qui va favoriser l’érosion des particules fines, d’où la 


récupération d’une eau de qualité mauvaise qui va nécessiter un post-traitement. 


v L'élévation en largeur à l’entrée : 


484,84 


484,80 


484,72 
484,68 


434,64 


Elevation 


434,60 


484,56 


484,52 





Chainage 


Figure 30: L'élévation du bassin en largeur à l'entrée après saturation 


On constate qu’il y a une irrégularité concernant la distribution de la boue en largeur, ce 


problème est dû au jet dispersif à partir des 3 entrées et leur mauvaise disposition. 





LE, A 


Mémoire du projet de fin d'études =/ N 


2. Scenario 2 : Digue d’épandage à pente 1%o sans chicanes 


v Géométrie : 


DTM View 3DView Results Verification Results 





Figure 31: géométrie de la digue à une pente 1%o 





l I | 
Tableau 12: résultats de récupération d’eau de la digue à une pente 1%o 


On constate que cette configuration du bassin n’est pas bénéfique vu que le taux de 
récupération ne dépasse pas 41%. Concernant la durée de la 1°° récupération, on constate 
qu’elle est bien inférieure à celle de la pente 0%, c’est normal vu l’existence d’une pente 


accélérant l’écoulement vers les points de collecte. 


V Le bassin après saturation : 


Mémoire du projet de fin d'études =/ W 


LEE, A 


ECOLE MOHAMMADIA D'INGÉNIEURS 


DTM View  3D View Results Verification Results 


Summary x 
183 000 
Date: 24 févr. 2024 
182 900 Raise Volume: 179 675 
Raise Tonnage: 228 621 
ne Cumulative Volume: 6 309 713 
Cumulative Tonnage: 8 109 415 
182 700 
182 600 
182 500 
182 400 
182 300 
182 200 
182 100 
182 000 
181 900 
8 8 8 8 8 8 
` Ñ ` Ñ Ñ Ñ 





Figure 32: l'élévation de la boue après saturation de la digue à une pente 1%o 
Date: 24 févr. 2024 
Raise Volume: 179 675 
Raise Tonnage: 228 621 
Cumulative Volume: 6 309 713 


Cumulative Tonnage: 8 109 415 








Tableau 13: taux de stockage de la digue à une pente 1 
On constate que cette configuration nous donne une augmentation du taux de stockage par 
rapport à la première configuration, et c’est dû à l’effet de la pente qui donne une énergie 

cinétique aux particules solides pour avancer vers l’avant et gagner plus de volume dans le 


bassin. 


Mémoire du projet de fin d’études =/ AN 
EE, A 


ECOLE MOHAMMADIA D'INGÉNIEURS 





Tableau 14: temps de séjour de l’eau pour une pente 1%o 


V Les formules utilisées : 


longueur du bassin (m) 


temps de séjour en (s) = ——— — — 
p J 2 vitesse moyenne d'écoulement (m/s) 


Le temps de séjour est un paramètre très intéressant, 1l nous donne une valeur moyenne 
concernant la durée qu’une goutte d’eau passe au sein du bassin. Tant que ce paramètre est 
grand, tant que l’eau récupéré à la sortie est de bonne qualité. Pour ce cas, on a trouvé 


11h27min, cette valeur est bonne ce qui nous mène à dire que l'eau récupéré est claire. 


v PL’'élévation du bassin en longueur après saturation : 


Elevation 





100 200 300 400 500 00 700 800 J300 1000 1100 1 200 
Chainage 


Figure 33: L'élévation du bassin en longueur après saturation à une pente 1%o 


On constate que également que la boue consolidée et drainée n’est pas bien répartie, en fait, la 
boue se décante sur tout le bassin mais 1l y a toujours ce problème de décantation auprès des 


entrées, ce qui explique encore la valeur médiocre du pourcentage de stockage. 


v L'élévation en largeur à l’entrée : 





Mémoire du projet de fin d’études =/ AN 


A 
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Figure 34: L'élévation du bassin en largeur à l'entrée après saturation 1%o 


On constate également qu’il y a une irrégularité concernant la distribution de la boue en 
largeur, ce problème est dû au jet dispersif à partir des 3 entrées et leur mauvaise disposition. 


Les trois sommets correspondent aux trois entrées fixées. 


v  L’élévation en largeur à la sortie : 


Elevation 





Chainage 


Figure 35: l'élévation en largeur à la sortie à une pente 1%o 


La cote de la boue à la sortie arrive à 475 NGM et le terrain naturel est à 470,81 NGM, donc 


on aura une élévation de 4,19 m de solide à la sortie. 





Mémoire du projet de fin d'études =, A 


pe) A 


ECOLE MOHAMMADIA D'INGÉNIEURS 


mm) Ce scenario est mieux que le premier, mais il n’est pas encore optimal, le 


pourcentage de stockage est encore faible ainsi que le pourcentage de récupération. 


3. Scenario 3 : Digue d’épandage à pente 1% sans chicanes 


v Géométrie : 


DTM View 3D View Results Verification Results 





Figure 36: géométrie de la digue à une pente 1% 





Figure 37: résultats de récupération d'eau de la digue à une pente 1% 


On constate également que cette configuration du bassin n’est pas bénéfique vu que le taux de 
récupération ne dépasse pas 52%. Concernant la durée de la 1°° récupération, on constate 
qu’elle est bien inférieure à celle de la pente 0%, c’est normal vu l’existence d’une pente 


accélérant l’écoulement vers les points de collecte. 


V Le bassin après saturation : 


A 


s . . . 2 s v 
Mémoire du projet de fin d’études À 
EE, 
ECOLE MOHAMMADIA D'INGÉNIEURS 
DTM View 30 View Results Verification Results 
Summary x 
Date: 16 juin 2025 
Raise Volume: 204 934 
Raise Tonnage: 258 678 
Cumulative Volume: 8 102 017 
Cumulative Tonnage: 10 297 287 


183100 
183000 
182900 
182700 
182600 
182500 
183400 
182300 
187200 
182100 
182000 
181900 
18188 

Ñ 





272 250 
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274 000 


Figure 38: l'élévation de la boue après saturation de la digue à une pente 1% 


Date: 16 juin 2025 
Raise Volume: 204 934 
Raise Tonnage: 258 678 

Cumulative Volume: 8 102 017 


Cumulative Tonnage: 10 297 287 


13 490 708,93 | 8 102 017,00 | 60,06 


Tableau 15: taux de stockage de la digue à une pente 1% 





On constate que cette configuration nous donne une augmentation du taux de stockage par 
rapport à la première configuration, et c’est dû à l’effet de la pente qui donne une énergie 
cinétique aux particules solides pour avancer vers l’avant et gagner plus de volume dans le 


bassin. 





Mémoire du projet de fin d'études =/ N 
LEE, 7 \ 


ECOLE MOHAMMADIA D'INGÉNIEURS 





Tableau 16: temps de séjour de l’eau pour une pente 1% 


Le temps de séjour a diminué vu l’existence de la pente. En fait, cette valeur est très faible, et 


on conclut qu'il y'aura une érosion forte sur tout le bassin et l'eau récupérée est turbide en 


particules solides fines. 


v  L’élévation du bassin en longueur après saturation : 
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Figure 39: L'élévation du bassin en longueur après saturation 1% 


On constate également que la boue consolidée et dramée n’est pas bien répartie, en fait, la 


boue se décante sur tout le bassin mais 1l y a toujours ce problème de décantation auprès des 


entrées, ce qui explique encore la valeur médiocre du pourcentage de stockage. 


v L'élévation en largeur à l’entrée : 





KA 


Mémoire du projet de fin d’études =/ AN 


435,00 


484,95 


434,90 
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Figure 40: L'élévation du bassin en largeur à l'entrée après saturation 1% 


On constate également qu’il y a une irrégularité concernant la distribution de la boue en 


largeur, ce problème est dû au jet dispersif à partir des 3 entrées et leur mauvaise disposition. 


Les trois sommets correspondent aux trois entrées fixées. 


Y D’élévation en largeur à la sortie : 


Elevation 





Chainage 


Figure 41: l'élévation en largeur à la sortie à une pente 1%o 


La cote de la boue à la sortie arrive à 467 NGM et le terrain naturel est à 460,17 NGM, donc 


on aura une élévation de 6.83 m de solide à la sortie. 





Ka 


Mémoire du projet de fin d'études ==, AN 
AY} IN 


m Ce scenario est mieux que les deux premiers, mais il n’est pas encore optimal, le 


pourcentage de stockage et de récupération sont encore faible. 


4. Scenario 4 : Digue d’épandage à pente 2% sans chicanes 


v Géométrie : 





Figure 42: géométrie de la digue à une pente 2% 





Tableau 17: résultats de récupération d ‘eau de la digue à une pente 2% 


On constate que cette configuration du bassin nous donne un taux de récupération intéressant. 
Concernant la durée de la 1™° récupération, on constate qu’elle est bien inférieure à celle de la 
pente 0%, c’est normal vu l’existence d’une pente accélérant l’écoulement vers les points de 


collecte. 


V Le bassin après saturation : 


Mémoire du projet de fin d'études =/ W 


pe) A 


ECOLE MOHAMMADIA D'INGÉNIEURS 


DTM View  3D View Results Verification Results 


Summary x 
Date: 27 mars 2030 
182900 Raise Volume: 203 54 
Raise Tonnage: 257 584 
182800 
Cumulative Volume: 14630 518 
182700 Cumulative Tonnage: 18 275 425 
182600 
182500 
182400 
182300 
182200 
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181900 
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Figure 43: l'élévation de la boue après saturation de la digue à une pente 2% 


Date: 27 mars 2030 
Raise Volume: 203 546 
Raise Tonnage: 257 584 

Cumulative Volume: 14 630 518 
Cumulative Tonnage: 18 275 425 





Tableau 18: taux de stockage de la digue à une pente 2% 


On constate que cette configuration nous donne une augmentation du taux de stockage par 
rapport à la première configuration, et c’est dû à l’effet de la pente qui donne une énergie 
cinétique aux particules solides pour avancer vers l’avant et gagner plus de volume dans le 


bassin. 





Mémoire du projet de fin d’études =/ AN 
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Tableau 19: temps de séjour de l’eau pour une pente 2% 


Le temps de séjour a diminué vu l’existence de la pente. La valeur du temps de séjour a 


encore diminuée donc l'eau récupérée à la sortie est de mauvaise qualité. 


V L’élévation du bassin en longueur après saturation : 
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Figure 44: L'élévation du bassin en longueur après saturation 2% 


On constate que cette configuration est très intéressante en terme de stockage, on arrive à 


stocker la boue sur les 4/5 du bassin (environ 80% du total du bassin). 


v L'élévation en largeur à l’entrée : 


Elevation 
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434,90 
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Figure 45: L'élévation du bassin en largeur à l'entrée après saturation2% 


On constate également qu’il y a une irrégularité concernant la distribution de la boue en 
largeur, ce problème est dû au jet dispersif à partir des 3 entrées et leur mauvaise disposition. 


Les trois sommets correspondent aux trois entrées fixées. 


Y  L'’élévation en largeur à la sortie : 
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Figure 46: l'élévation en largeur à la sortie à une pente 2% 


La cote de la boue à la sortie arrive à 475 NGM et le terrain naturel est à 448,35 NGM, donc 
on aura une élévation de 26,65 m de solide à la sortie. Ce qui montre la répartition plus ou 


moins bonne de la boue sur tout le volume du bassin. 


mm} Ce scenario est mieux que les premiers, on peut dire qu’il est optimal, le pourcentage 


de stockage et de récupération sont bons et le temps de séjour est moyen. 


+Æ on constate que pour les scenarios précédents, on avait le problème de la 
pente raide en longueur qui engendre une récupération de l’eau d’une qualité 
mauvaise, d’où on a pensé à introduire des chicanes qui vont jouer le rôle 
d’un obstacle afin de limiter l’effet de la pente et augmenter le temps de 


séjour pour avoir une eau claire à la sortie. 


5. Description des chicanes utilisées : 





Ka 


LEE, A 
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hi 6 PAL] 0,30 236,48 3 


Tableau 20: description des chicanes utilisées 





Figure 47: vue en 3D d'une chicane 


6. Scenario 5 : Digue d’épandage à pente 1%0 avec introduction des 


chicanes 


vV Géométrie : 


DTM View  3D View Results Verification Results 
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Figure 48: géométrie de la digue à une pente 1%o avec chicanes 





Figure 49: résultats de récupération d’eau de la digue à une pente 1%o avec chicanes 


On constate que cette configuration du bassin nous donne un taux de récupération meilleur 
que celui trouvé pour la même pente et sans chicanes, mais il est encore médiocre. L’effet des 
chicanes se manifeste en temps de la 1™° récupération qui a augmenté de 11 jours par rapport 


à la même pente et sans introduction des chicanes. 


V Le bassin après saturation : 


DTM View  3D View Results Verification Results 
Summary x 
Date: 29 sept. 2023 
Raise Volume: 161 141 
Raise Tonnage: 206 238 
Cumulative Volume: 5 778 278 
Cumulative Tonnage: 7 433 242 
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Figure 50: l'élévation de la boue après saturation de la digue à une pente 1%oavec chicanes 


Date: 29 sept. 2023 
Raise Volume: 161 141 
Raise Tonnage: 206 238 

Cumulative Volume: 5 778 278 


Cumulative Tonnage: 7 433 242 
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Tableau 21: taux de stockage de la digue à une pente 1 %o avec chicanes 


On constate que cette configuration nous donne un taux de stockage inférieur à celui trouvé 
pour la même pente sans chicanes 66%, ça s’explique par l’effet d’obstacle joué par les 
chicanes, ces derniers diminuent la dispersion des particules solides, d’où la perte en stockage 


dans les 600 m derniers. 


Ÿ Longueur approximative parcourue: 


DTM View  3D View Results Verification Results 





Summary 
183 000 X 
Date: 16 mars 2021 
182 900 Raise Volume: 63 572 
Raise Tonnage: 86 072 

182 800 

Cumulative Volume: 2 497 206 
182 700 Cumulative Tonnage: 3 194 398 
182 600 
182 500 
182 400 
182 300 


w229 Longueur totale: 


w20 Lr= 5*l+ 4*L= 3408m 


182 000 
181 900 

8 8 8 8 8 8 
RS: R R R g Š 








Tableau 22: temps de séjour de l’eau pour une pente 1 %o avec chicanes 


Le temps de séjour a augmenté par rapport à la configuration 1%o sans chicanes, et c’est dû à 
P introduction des chicanes qui augmentent le temps qu’une goutte d’eau passera dans le 


bassin afin d’avoir une auto-purification et récupérer une eau claire aux points de collecte. 


Mémoire du projet de fin d’études =/ AN 


LE} | A 


ECOLE MOHAMMADIA D'INGÉNIEURS 


v  L’élévation du bassin en longueur après saturation : 


DTM View  3D View Results Verification Results À 
Summary Xx 
Date: 29 sept. 2023 
182900 Raise Volume: 161 141 
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Figure 52: L'élévation du bassin en longueur après saturation 1%o 


> suivant la ligne contenant la premiére entrée coupe A-A’ : 
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Figure 53: l'élévation en longueur du bassin après saturation à une pente de I %o avec chicanes 


On constate que cette configuration est très intéressante en terme de diminution des vitesses 
par les chicanes, donc on a résolu le problème de la pente raide de 1%, mais on a diminué 
également l’aspect dispersif de solide, ce qui explique la partie presque non exploitable par la 


boue après la chicane sur le graphe. 
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YV L’élévation en largeur à l’entrée coupe B-B’ : 
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Figure 54: L'élévation du bassin en largeur à l'entrée après saturation 1%o 


On constate également qu’il y a une irrégularité concernant la distribution de la boue en 
largeur à l’entrée, ce problème est dû au Jet dispersif à partir des 3 entrées et leur mauvaise 


disposition. Les trois sommets correspondent aux trois entrées fixées. 


Y  L'’élévation en largeur à la sortie coupe C-C’ : 


Elevation 





Chainage 


Figure 55: l'élévation en largeur à la sortie à une pente 1 %o 


Cette coupe nous explique pourquoi le taux de récupération pour cette configuration est 
faible, on constate qu’un volume non négligeable va être stagnant dans la partie gauche du 


bassin loin des points de récupération (voir les points de collecte sur les plans ci Joint). 
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m Ce scenario est à éviter. 


7. Scenario 6: Digue ď’épandage à pente 1% avec introduction des 


chicanes 
v Géométrie : 


DM view Den Ress Verification Results 





Figure 56: géométrie de la digue à une pente 1% 





Tableau 23: résultats de récupération d'eau de la digue à une pente 1% avec chicanes 


On constate que cette configuration du bassin nous donne un taux de récupération meilleur 
que celui trouvé pour la même pente et sans chicanes. On constate également la diminution du 


temps de la 1° récupération. 


V Le bassin après saturation : 
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Figure 57: l'élévation de la boue après saturation de la digue à une pente 1% avec chicanes 


Date: 14 juin 2025 
Raise Volume: 179 849 
Raise Tonnage: 228 112 
Cumulative Volume: 8 066 950 
Cumulative Tonnage: 10 286 635 





Tableau 24: taux de stockage de la digue à une pente 1% avec chicanes 


On constate que cette configuration nous donne un taux de stockage proche à celui trouvé 


pour la même pente sans chicanes 60%. 





Tableau 25: temps de séjour de l'eau pour une pente 1% avec chicanes 
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Le temps de séjour a augmenté par rapport à la configuration 1% sans chicanes, et c’est du à 


l'introduction des chicanes qui augmentent le temps qu’une goutte d’eau passera dans le 


bassin afin d’avoir une auto- 


purification et récupérer une eau claire aux points de collecte. 


V L’élévation du bassin en longueur après saturation : 
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Figure 58: L'élévation du bassin en longueur après saturation 1% avec chicanes 


Le graphe nous montre la cause derrière le faible pourcentage de stockage pour cette 


configuration, en fait le volume derrière la dernière chicane n’est pas exploité. 


Y L’élévation en largeur à l’entrée: 
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Figure 59: L'élévation du bassin en largeur à l'entrée après saturation 1% avec chicanes 
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On constate qu’on a toujours ce problème de la répartition non uniforme en largeur de la 
boue à l’entrée du bassin, ce problème va être résolu après quand on trouvera la meilleur 


configuration en terme de stockage et récupération de l’eau. 


Y  L'’élévation en largeur à la sortie : 
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Figure 60: l'élévation en largeur à la sortie à une pente 1% avec chicanes 


Le problème de la répartition non uniforme de la boue à la sortie du bassin est résolu par la 
multiplication des points de collecte, en fait il y a 3 points de collecte sur le côté sud du bassin 


espacés de 100 m et le quatrième existe sur le côté nervure. 


8. Scenario 7 : Digue d’épandage à pente 2% avec introduction des 


chicanes 


V Géométrie : 
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Figure 61: géométrie de la digue à une pente 2% avec chicanes 





Figure 62: résultats de récupération d’eau de la digue à une pente 2% avec chicanes 


Cette configuration nous donne un pourcentage très élevé concernant la récupération de l’eau 
et la durée de saturation est la meilleure parmi les scenarios étudiés, la durée de la 17° 


récupération est bonne en comparaison avec les autres scenarios. 


Mémoire du projet de fin d’études =/ AN 


LE} A 


ECOLE MOHAMMADIA D'INGÉNIEURS 


V Le bassin après saturation : 


DTM View  3D View Results Verification Results 
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181900 


Summary x 
Date: 10 sept, 2030 
Raise Volume: 204 169 
Raise Tonnage: 258 248 
Cumulative Volume: 15 223 467 
Cumulative Tonnage: 19 041 262 
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Figure 63: l'élévation de la boue après saturation de la digue à une pente 2% avec chicanes 


Date: 10 sept. 2030 
Raise Volume: 204 169 
Raise Tonnage: 258 248 

Cumulative Volume: 15 223 467 


Cumulative Tonnage: 19 041 262 





19 211 431,80 15 223 467,00 





Tableau 26: taux de stockage de la digue à une pente 2% avec chicanes 


On constate que le pourcentage de stockage est intéressant en comparaison avec les scenarios 


précédents. 
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Tableau 27: temps de séjour de l’eau pour une pente 2% avec chicanes 


Y L'élévation du bassin en longueur après saturation : 
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Figure 64: L'élévation du bassin en longueur après saturation 2% avec chicanes 


On constate que la boue prend une répartition uniforme en longueur et que les chicanes sont 
toutes couvertes de la boue ce qui montre l’aspect dispersif fort de la boue renforcé par la 


pente. 


v L'élévation en largeur à l’entrée : 
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Figure 65: L'élévation du bassin en largeur à l'entrée après saturation 2% avec chicanes 


Mémoire du projet de fin d’études —/ \ \ 


ECOLE MOHAMMADIA D'INGEÉNIEURS 


On constate que le problème de la répartition non uniforme de la boue à l’entrée existe 


toujours. 


Y L’élévation en largeur à la sortie : 
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Figure 66: l'élévation en largeur à la sortie à une pente 2% avec chicanes 


On constate que la répartition de la boue est uniforme à la sortie du bassin et l’élévation 


importante de la boue qui arrive à 27m de hauteur par rapport au terrain naturel. 


9. Conclusion : 
À partir des statistiques calculées, on conclut que le meilleur scenario est celui de la 
configuration du bassin avec une pente de 2% avec introduction des 4 chicanes, on constate 
également que la répartition de la boue en longueur est bonne mais la répartition de la boue à 
l’entrée du bassin en largeur n’est pas uniforme, c’est pourquoi on a pensé d’ajouter 2 entrées 
au droit des sections non saturées après la durée de saturation prévue pour cette configuration 


(voir le graphe suivant). 
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Figure 67: positions d'ajout des 2 entrées côté Nord 


10. Amélioration du modèle optimal : 


v Résultats de l’ajout des 2 entrées : 
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Figure 68: résultats de l'ajout des 2 entrées Nord sur l'élévation du bassin 


Date: 11 sept. 2030 


Raise Volume: 2 984 


Raise Tonnage: 4 168 


Cumulative Volume: 5 184 





Cumulative Tonnage: 7 242 


On constate que l’ajout des 2 entrées n’a pas eu un effet important concernant le volume de 
solide stockés vu que la hauteur à l’entrée où il y’avait le problème était de 0.5 m et sur une 


distance de 250m, après cette distance en longueur on trouve une pente bien uniforme. 
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v Coupe en largeur le 11 septembre 2030 : 


Figure 69: résultats de l'ajout des 2 entrées sur l'élévation du bassin en largeur à l'entrée 


Elevation 


Vu qu’on ne gagne pas un volume intéressant en terme de stockage ou bien en terme de 
récupération d’eau, l’ajout des 2 entrées au côté nord du bassin est à éviter, et on va 
essayer l’ajout des entrées aux cotés EST et OUEST à une distance de 591m à partir du 


côté Nord du bassin c’est-à-dire à m1-distance en longueur. 


Ÿ Les coordonnées des points de déposition 


Tableau 28: coordonnées des 2 entrées Est et Ouest ajoutées 





273311,4698 182442,8626 
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DTM View  3D View Results Verification Results 
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Figure 70: vue d'emplacement des 2 entrées 


v Présentation des résultats : 





Tableau 29: résultats d'ajout des 2 entrées Est et Ouest 


On constate que l’ajout des entrées nous a permet de gagner en terme d’exploitation de la 
digue d’épandage environ 2 ans de plus, et le taux de récupération qui est amélioré de 88,8% 


à 92,56%. 
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V Le bassin après saturation : 
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Figure 71: lélévation de la boue après saturation de la digue et ajout des 2 entrées 


Date: 05 août 2032 
Raise Volume: 44 209 
Raise Tonnage: 59.52? 


Cumulative Volume: 2 412 428 


Cumulative Tonnage: 3 180 215 





Tableau 30: taux de stockage de la digue à une pente 2% 


On constate que le pourcentage de stockage pour la quantité ajoutée par les 2 entrées est 


intéressant. Maintenant on calcule le pourcentage de stockage total : 
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Tableau 31: taux de stockage de la digue à une pente 2% après ajout des 2 entrées 


On constate que le pourcentage total de stockage a augmenté de 76,16% avant l’ajout des 2 


entrées à 88,71% après l’ajout des 2 entrées. 





Tableau 32: temps de séjour de l’eau pour une pente 2% après ajout des 2 entrées 


Le temps de séjour est faible vu la distance courte parcouru entre les 2 entrées et les points de 


collecte et la pente assez raide construite après déposition (pente de 2%). 
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Figure 72: vue des coupes réalisées 
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V L’élévation du bassin en longueur à partir des deux entrées coupe A-A’ : 
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Figure 73: coupe A-A' 


On constate que la boue prend une répartition uniforme en longueur et qu’on atteint une cote 


NGM intéressante à la sortie de 480NGM. 


Y L’élévation en largeur à la section des 2 entrées coupe B-B’ : 
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Figure 74: coupe B-B' 


On constate que la répartition non uniforme de la boue en largeur existe toujours mais 


puisqu’il y a une pente de 1% en longueur qui favorise l’écoulement vers la sortie, on ne 


craint rien. 


Ka 
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Y L’élévation en largeur à la sortie coupe C-C’ : 
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Figure 75: coupe C-C' 


On constate que la répartition de la boue est uniforme à la sortie du bassin et l'élévation 


importante de la boue qui arrive à 33m de hauteur par rapport au terrain naturel. 


Remarque : Coupe D-D’ 
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Figure 76: coupe D-D' 


On constate qu’on va perdre un volume d’eau qui va être emprisonné dans l’espace indiqué 


sur le graphe, mais il apparait qu’il est négligeable. 


11. Recommandations et planning dďd’utilisation du bassin : 


Pour conclure, le meilleur scenario est de construire un bassin avec les instructions suivantes : 





Ka 
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e La surface de base à une pente de 2%. 

e Introduction des 4 chicanes. 

e Respecter la distribution des 3 entrées au côté Nord du bassin. 

e Ajouter 2 entrées à m1-distance en longueur après la saturation. 

e Pour avoir les résultats présentés et l’efficacité en terme de rendement de stockage et 


de récupération pour cette configuration, 1l faut respecter le planning des quantités à 


débiter mentionné en annexe 1.]. 
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Figure 77: 5 entrées de boue proposées 


On obtient le planning en se basant sur les tonnages déposés dans chaque étape, fournis par 
RIFT TD, après en divisant par la durée de déposition, on trouve le débit solide de chaque 
entrée, puis en multipliant par le rapport débit solide/débit boue, on trouve le débit à desservir 


à partir de chaque entrée sur toute la durée de saturation. 
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Chapitre V : Stabilité de la digue et dimensionnement 
des chicanes 
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I. Stabilité de la digue: 

La zone d’étude constitue une ancienne zone d’exploitation des phosphates et 
renferme par conséquent des zones en tranchées généralement remblayées par le stérile du 
phosphate. A ceci, s’ajoute la présence d’une série de cavaliers. Une telle configuration a 
chamboulé la topographie et la géomorphologie de la future digue. Ainsi, le site d’étude 
correspond à une zone généralement en remblai avec une hauteur variable. 

Les calculs de stabilité ont été menés en considérant : 


e La géométrie de la fondation ; 
e Les caractéristiques mécaniques des différents matériaux constituant la digue ; 


e Les caractéristiques mécaniques de la fondation au droit de la digue. 


1. Paramètres de modélisation : 


Les paramètres suivants sont retenus pour l’ensemble des modélisations, ils sont 


obtenus à partir des différents essais géotechniques effectués par le laboratoire LPEE. 


Re 





Tableau 33: paramètres géotechniques utilisées dans la simulation 


2. Profil de calcul : 

L’étude de stabilité de la digue a été effectuée sur deux profils types (Les cas 
critiques en termes des hauteurs et des quantités remblai), le 1° profil est situé sur le côté 
Nord et constitué notamment de cavaliers existants avec une petite quantité de remblai en 
amont et en aval, tandis que le 2ème est situé sur le côté Sud constitué totalement en remblaï 


(Voir les figures 78 et 79). Ces profils sont symétriques avec une pente de 2.6H/I V. 
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3. Méthodologie de calcul : 
Les calculs de stabilité ont été effectués à l’aide des logiciels SEEP/W et SLOPE/W, 


faisant partie de la suite GeoStudio (version 2012) : 
V SEEP/W permet de calculer la ligne piézométrique et les pressions interstitielles pour 
les différents cas considérés, ces résultats sont utilisés comme données d’entrée dans 


le calcul de stabilité ; 


V SLOPE/W permet de calculer la stabilité des pentes suivant les méthodes analytiques 
citées antérieurement, comme il offre la possibilité de définir la région des centres de 
cercles de glissement et des rayons limites, ou bien d’effectuer la recherche 


automatique de la surface de rupture critique. 


La ligne piézométrique est calculée par éléments finis sur SEEP/W, un écoulement en régime 
permanent est considéré pour le remplissage par la boue accompagnée d’une couche 
supérieure d’eau, avec une perméabilité isotrope et homogène dans le corps de la digue et la 


fondation. 


Puisque l’écoulement transitoire d’un matériau n’est pas possible à modéliser (remplissage 
par la boue avec une couche supérieure d’eau), le calcul a été effectué en régime permanent 
pour plusieurs niveaux de remplissage de la boue. 

La stabilité de la digue est calculée par la méthode des tranches, basée sur la théorie de 
l’équilibre limite. La méthode consiste à examiner l’équilibre des forces actives agissant le 
long de la surface de rupture potentielle et les forces de résistance agissant le long de la même 
surface. 

Le calcul de stabilité est ramené à un problème plan en négligeant les interactions. Ce choix 
est justifié par la géométrie relativement simple de la digue. Ce choix est sécuritaire par 
rapport à un calcul 3D puisqu'il considère que la rupture se produit le long de la digue, et 
néglige l’interaction entre les sections. 

Dans le cadre de ce projet, la méthode de Morgenstern-Price est utilisée, et l’algorithme de 
recherche de surface critique est conduit en deux étapes : Recherche de surface critique 
circulaire par un processus itératif, puis l’optimisation de cette dernière pour obtenir une 
surface critique de forme quelconque. 


Le critère de rupture de tous les sols est Mohr-Coulomb. 
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4. Situations envisagées : 


1. Fin de construction : 
Le corps de la digue est considéré constitué de matériau assez étanche. En 


conséquence, et sous l’effet de la remontée des remblais, des pressions interstitielles sont 
générées dans le corps de la digue. 

Ces pressions interstitielles réduisent la stabilité de l’ouvrage à court terme, avant de se 
dissiper lentement par consolidation. La prise en compte de ce phénomène conduit à aborder 


le cas de charge de fin de construction. 


2. Retenue vide : 
Le cas d’une retenue vide est traité pareillement au cas de fin de construction, avec 


une modification des paramètres géotechniques du matériau de la digue, en raison de la 


dissipation lente des pressions interstitielles par consolidation. 


3. Retenue remplie par la boue à hcritique : 
Le niveau de la boue est traité en plusieurs cas afin de trouver la hauteur critique vis- 


à-vis la stabilité de la digue. Une couche infinitésimale d’eau est supposée existante en haut 
de la boue afin de modéliser la ligne piézométrique. Le calcul des conditions hydrauliques 


dans la digue est réalisé en régime permanent. 


4. Cas sismique : 
D’après RPS 2000 version 2011, la commune de Ben Guérir relevant de la province 


El Kelaâ Sraghna est située dans la zone 1 en accélération et en vitesse, ce qui correspond à 
une accélération de 0.07g, et une vitesse de 0.07 cm/s. La probabilité de dépassement de ces 
valeurs est égale à 10% en 50 ans, soit une période de retour de 475 ans. Par analogie avec le 
calcul des actions sur les bâtiments proposés par le RPS. L’accélération est multipliée par des 
facteurs dépendant du site, et de la classe de l’ouvrage. Les coefficients suivants sont 


proposés pour la digue de stockage des boues de Benguerir : 


v Coefficient de site : S1 = 1,00 
v Coefficient de priorité : Classe III I = 1,30 


Il en résulte une accélération de calcul = 0.091. 

L’accélération maximale considérée dans le calcul de stabilité est : Gmax= 0.1g. 

L’action sismique est prise en compte par la méthode pseudo-statique, elle consiste à 
schématiser l’action dynamique du séisme par un chargement statique : l’effort appliqué sur 


un barrage est égal à une fraction de la force de pesanteur. La force sismique pseudo-statique 


est calculée par la formule : E 
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Avec a coefficient pseudo-statique égal à : 


e On =2/3 = 0.67 pour la composante horizontale ; 


e av —1/3* on = 0.22 pour la composante verticale ; 
Le niveau d’eau, les sous-pressions, et les paramètres géotechniques sont analogues à ceux de 
la retenue vide. 
Les situations sismiques sont considérées pour les deux cas suivants : 


e Retenue vide : 
e Retenue remplie par la boue à la hauteur critique. 


5. Résultats de la modélisation : 


1. Profil dimensionnant type 1 (Nord) : Haim=13m ; Hav=3m 


La digue utilisée pour le côté Nord et partie centrale du bassin est constituée des 





cavaliers existants sur terrain, avec ajout de petites quantités de remblai de l’extérieur du 
bassin, et vu la qualité des matériaux de ces cavaliers, on n’a pas besoin d’ajouter un drain ou 


une géomembrane, le profil dimensionnant est le suivant : 


Fondation 





Figure 78: profil dimensionnant Nord et partie centrale du bassin 


Le tableau suivant résume les facteurs de sécurité critiques obtenus : 








Mémoire du projet de fin d'études =/ AN 
EE, 


A 


Tableau 34: résultats de stabilité du profil type 1 
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On constate que tous les cas étudiés sont vérifiés en terme de stabilité, les résultats détaillés 
sont joints dans l’annexe 2.1. On peut également augmenter le coefficient de sécurité par 
l’ajout d’une butée à l’aval de la digue afin de la stabiliser de plus, ou bien ajouter une 
géomembrane qui joue le rôle d’un obstacle devant l’eau ainsi elle diminue les pressions 
interstitielles qui se développe à l’intérieur de la digue, on peut également ajouter un drain 


afin d’annuler rapidement les pressions interstitielles développés. 


2. Profil dimensionnant type 2 (Sud) Ham=13m ; Hav=13,85m 


La digue utilisée pour le côté Sud du bassin est constitué de remblais transportés à 
partir des cavaliers de bonne qualité les plus proches (voir annexe). En fait, ce profil est 
dimensionné en ajoutant une (géomembrane +céotextile) à l’amont de la digue et un drain à 


son aval. Cet ajout est pour garantir une stabilité du profil Sud. 


Le profil inclut la stabilité de la tranchée qui collecte l’eau à l’aval de la digue et le dirige vers 


le bassin de récupération des eaux. Le profil dimensionnant est le suivant : 


du o A 





Figure 79: profil dimensionnant Sud 


Le tableau suivant résume les facteurs de sécurité critiques obtenus : 


profil sud 


Cas sismique 





LEE, A 
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Tableau 35: résultats de stabilité du profil type2 


On constate que tous les cas étudiés sont vérifiés en terme de stabilité, les résultats détaillés 


sont joints dans l’annexe 2.2. 


II. Dimensionnement des chicanes: 
Les chicanes introduites pour le cas optimal du bassin (pente 2% avec chicanes) sont 


considérées comme des murs de soutènement qui résistent à la poussée de la boue cumulée. 


1. Predimensionnement: 


D'après SETRA (service d'études sur les transports, les routes et leurs aménagements): 
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Figure 80: predimensionnement des murs de soutènement selon SETRA 


Hauteur de prédimensionnment de chicane: H= 7m 


e Epaisseur du voile en amont : 


Cv —(,3 m 
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0, 15<= Cv 


e Base de la semelle : 
B=3,5 m 
0,4H<=B<=0,6H 

e Patin 


P= 0,9 m 


B /5<=P<=B /2 
e Epaisseur du voile : 
Evoile = 0,6 m 
e Epaisseur de la semelle : 


Csemelle = 0,6 M = voile 
e Fruit: 


Fruit=5,5 % >= 2% 


Donc on résume notre prédimensionnement : 


Leh] 
ran 


NW 


KR 











Figure 81: Géométrie de prédimensionnement 


2. Inventaire des charges: 


e  Poussée et butée des sols : 


Pour la force de poussée : 





Mémoire du projet de fin d’études —/ \ AN 


ECOLE MOHAMMADIA D'INGÉNIEURS 


x = emf ©) 
4 2 


fr: 


Coefficient de poussée: 
La force de poussée résultante est : Fa = =KayH - 


Pour la butée: 


nE | K, = tan (3+ ai 
Coefficient de butée: F 4 2 


La force de butée résultante : Fp = ZK pyH 


Soit les notations suivantes: 

N1 : Poids des terres en amont 

N2 : Poids des terres en aval 

N3 : Poids du mur 

N4 : Poids de la semelle 

Fa : Résultante de la force de poussée 


Fp : Résultante de la force de butée 


a ž Fa 


N4 


Figure 82: notations des forces appliquées 


3. Modélisation sur CYPE: 
En utilisant les données géotechniques des sols déjà citées dans la partie de stabilité de la 
digue dans le tableau 34, on modélise sur CYPE notre mur de soutènement avec les forces des 


terres en suivant la méthode de Rankine. 
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30cm 
0.00 m (cm) a E 0.00 m 





-5,40 m -5,40 m 
-7.00 m 
-7.60 m -7.60 m 


papi mha obus — m (em) 
Figure 83: Modélisation des forces des terres 
4. Résultats de dimensionnement en CYPE: 
Après introduction de la géométrie et les caractéristiques du sol, j'ai trouvé que toutes les 
conditions de stabilité sont vérifiées sauf la condition de non renversement, alors pour 


remédier à ce problème de conception, j'ai élargis le patin de 90cm à 150 cm afin d'avoir un 


moment résistant au renversement par rapport au point B (m1-distance de la semelle). 


Le profil vérifiant les conditions de stabilité est le suivant : 
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0.00 m (cm) -a 0.00 m 





. . - 
’. " ' . 
t g ' 
E C A 


-5.40 m -5.40 m 
-7.00 m 
-T.60 m -7.60 m 


| bi a — {cm} 
Figure 84: Géométrie finale du mur de soutènement 
1. Stabilité au glissement : 


Pour vérifier la condition de stabilité vis-à-vis le glissement, il faut que : 


E,*tano' 
PTT 0 
Fy 


Avec @" = T (angle de frottement du sol). 


F,: la somme des forces verticales stabilisantes. 
Fg: la somme des forces horizontales déstabilisantes. 


Le cercle de glissement critique trouvé est le suivant : 
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21.6 m ———————— 


1215 i E a. #45 10 E 
11.00 m i i i i , l 





Terrain Naturel 
0.00 m EX can 0.00 m 


-5.40 m -5.40 m 


-7.00 m A 
-7.60 m RE _ -7.60 m 


Figure 85: Cercle de glissement critique 


Avec un coefficient de sécurité =2,35 conforme : 


Centre du cercle de glissement le 
plus défavorable | 
Coeficient de sécurité: 2.355 (CONFORME) 


Distance 3 l'arase du mur X: 0.49 m 
Cote Y- 0.78 m 
Rayon: 8.76 m 





2. Stabilité au renversement: 


Pour avoir la stabilité vis-à-vis le renversement, 1l faut que: 


M 
V/A > 1,5 
My ,A 


Avec M,;4: somme des moments des forces verticales par rapport à A 
MA: somme des moments des forces horizontales par rapport à À 


La note de calcul de CYPE nous donne: 


ZA 


A 
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Vérification Valeurs État 
Vérification de la stabilité: 
- Coefficient de sécurité au renversement: Minimum: 1.8 

Calculé: 2.02 Conforme 
- Coefficient de sécurité au glissement: Minimum: 1.5 

Calculé: 2.35 Conforme 
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3. Stabilité au poinçonnement: 


La condition à respecter est : 


ref = qadm 
Avec qref: contrainte conventionelle de référence (dépend de l'excentricité et 
du chargement de la géometrie) en kPa 
adm: contrainte admissible du sol en kPa 


CYPE nous donne la vérification de stabilité vis-à-vis le poinçonnement : 





Contraintes appliquées au sol: 


-Contrainte moyenne: Maximum: 0.2 MPa 
Calculé: 0.1278 MPa (Conforme 
-Contrainte maximale: Maximum: 0.25 MPa 


Calculé: 0.2468 MPa (Conforme 

















AVEC qadm= 0,25 MPa et Grer= 0,246 MPa 
4, Détermination des armatures: 
Selon le BAEL 99 


e Calcul de voile: 


Le voile est calculé comme une console encastrée sur la semelle et soumise à une charge 
trapézoïdale, les armatures et leurs arrêts se déterminent à partir des diagrammes du moment 


de flexion composée et de l’effort tranchant. 


Le dimensionnement se fait à l’état limite ultime en suivant l'organigramme mentionné en 


annexe 2.3 
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CYPE nous donne le ferraillage suivant: 


AVANT h ARRIÈRE 


S8FAHATĒe=20 
p 





JAP4HA16e=15 





Figure 86: Ferraillage du mur de soutènement 


5. Etanchéité de la chicane: 
Vu que le béton est très fragile en contact avec l'eau, et pour ne pas avoir un développement 
des pressions interstitielles à l'intérieur du mur qui peuvent s'infiltrer à l'intérieur du corps de 
la digue et toucher ainsi à sa stabilité, on prévoit couvrir les chicanes avec une 


(géomembrane+séotextile) afin d'imperméabiliser les chicanes contre l'infiltration des eaux. 
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Conclusion 


En guise de conclusion, l'étude du comportement de la boue à l'intérieur du bassin 
pour plusieurs scénarios proposés a montré que le cas optimal, en terme de stockage de solide 
et récupération d'eau claire, est de construire un bassin de pente 2% avec introduction des 


chicanes de hauteur 5 à 7m. 


Le dimensionnement des chicanes s'est fait en utilisant CYPE, et concernant la 
stabilité du corps de la digue pour le cas optimal, on doit disposer des quantités de remblais 
nécessaires pour la construction de la digue ainsi que ses profils dimensionnants. La 
vérification de la stabilité faite dans ce rapport concerne le premier scénario (bassin sans 


pente et sans chicanes) choisi par le bureau d'étude. 
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Annexe : 


I. Annexe 1: 


a. Repérage de la digue en google maps: 








b. Route d'accès à la digue d'épandage 
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c. Vue en plan de la digue d'épandage (géomembrane en gris) 


BASSIN DE FEPAISE 


~ CLARINETTE CONDUITE 
NMETALLIQUE CEDO mm 


~ CLARINETTE CONDUITE 
METALLIQUE 600 mm 


TAPIS DRAINANT 


~ GLARINETTE CONQUITE 
METALLIQUE CEDO mmi 


wpy'9SS =1nə8127 


TAPIS DRAINANT 


= TAPIS DRAINANT 


PERIMETRAL — GLARIMETTE CONDUME 
3 5" METALLIQUE E0 mm 


B — ~ mmm 


— 





~ UMME FOUILLES 


T GEOMEMBRANE ET GEOTEXTILE 


d. Détail de pose de la géomembrane 


DETAIL POSE SACS DE SABLES 


Géomembrane 


























~ Géotextile 





























2 Sacs de sable 
Tige en fer HA10 Espacement latéral=2.5m 
pour fixation 


de la géomenbrane 
Espacement latéral = 2.5m FR 
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ECOLE MOHAN 


Coupe longitudinale des cheminées drainantes 


€. 


"F= 


Flague mëtalligue 2 


mm: 


Ep= 





Plaque métallique N°1: 
Ep= mm 


3 


Exaler 


Aai 





wadi 
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f. Le nivellement de la zone d'implantation de la digue 





2. 
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h. Coupe lithologique 


Log stratigraphique Formations 


‘evahertrembhla compcse de blocs de colcaire, marn 


nosphèle non espiolte et autres elements decapes 


Marne jaunâtre phosphatrs 


Phos hate 9g eseis A «x Île M ONSatre 


Calcaire à sile» 


Phosphate gréseux à couleur grisétre 





1. Résultats d'essais de perméabilité Lugeon et Lefranc 
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j. Planning des débits des 5 entrées sur la durée de vie du bassin 


débit boue tonnes/heure 


débit entrée 1 | débit entrée 2 | débit entrée 3 | débit entrée 4 | débit entrée 5 
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Emi S 


13 
09/04/2032 14 


22/04/2032 


06/05/2032 14 
19/05/2032 13 
01/06/2032 13 


14/06/2032 13 


27/06/2032 
10/07/2032 


23/07/2032 13 


05/08/2032 
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II. ANNEXE 2: 


1. Résultats de stabilité du profil dimensionnant Nord: 


a. Cas Fin de construction : 


La ligne piézométrique 


Charge hydraulique totale 


Bm<0-2m 
O2-4m 


altitude (m) 





0 10 20 30 40 50 60 70 80 
distance (m) 


Facteur de sécurité talus amont 


10 


altitude (m) 


-10 





distance (m) 
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altitude (m) 


altitude (m) 
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Facteur de sécurité talus aval 





z0 
10 
0 
-10 
-20 
0 10 20 30 40 50 50 TŪ 20 
distance (m) 
b. Cas Retenue vide : 
Facteur de sécurité talus amont 
a 
10 
[A] 
-100 
-31 





A] LL] 21 sil 2] 1 aa] F0 i 


distance (m) 
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Facteur de sécurité talus aval 


altitude (m) 





distance (m) 


c. Cas Retenue remplie de boue à huritique =12m: 
La ligne piézométrique 


Charge hydraulique totale 


B = 10-10,2 m 
o 10,2-10,4 m 
o 10,4-10,6 m 
E 10,6 - 10,8 m 
m 10,8-11m 

m 11-11,2 m 

o 11,2-11,4 m 


oO 11,4-11,6 m 
o 11,6-11,8 m 
o 118-12m 
B 12-12,2 m 
E 212,2 m 






altitude (m) 


distance (m) 
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Facteur de sécurité talus amont 


20 


altitude (m) 


-20 
(ij 10 20 30 40 5o 


distance (m) 


Facteur de sécurité talus aval 


20 


10 


altitude (m) 


-10 


-20 
a 10 20 30 40 50 


distance (m) 


AS 
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d. Cas Retenue boue à hcritique =3m: 
La ligne piézométrique 


Charge hydraulique totale 


B <3,34-3,345 m 
ET 3,345 - 3,35 m 
M 3,35- 3,355 m 
E 3,355 - 3,36 m 
oO 3,36 -3,365 m 
E 3,365 -3,37 m 
E 337-3375 m 
T 3,375 -3,38 m 
A 3,36 -3,365 m 
o 3,385 -3,39 m 
o 3,39 - 3,395 m 
o 3,395 - 3,4 m 
O 3,4-3,405 m 
O 3,405 -3,41m 
oO 3,41 -3,415 m 
o 3,415 -3,42m 
E 3,42- 3,425 m 
B 3,425 -343 m 
B 2343m 


altitude (m) 





distance (m) 


Facteur de sécurité talus amont 


10 


altitude (m) 


-10 





0 10 20 30 40 50 g0 rO 20 


distance (m) 
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Facteur de sécurité talus aval 


altitude (m) 


30 40 


distance (m) 


e. Cas sismique retenue vide: 


Facteur de sécurité talus amont 


a 10 


altitude (m) 


30 40 50 
distance (m) 


0 10 20 
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AS 


TŪ 20 
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Facteur de sécurité talus aval 


altitude (m) 





0 10 20 30 40 50 60 F0 20 


distance (m) 


f. Cas sismique retenue de boue à hcritique =3M : 


Facteur de sécurité talus amont 


10 


altitude (m) 


-10 





0 10 20 30 40 50 a0 ra 20 


distance (m) 
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altitude (m) 





altitude (m) 
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Facteur de sécurité talus aval 


20 30 40 50 


distance (m) 


2. Résultats de stabilité du profil dimensionnant Sud: 


a. Cas Fin de construction : 


La ligne piézométrique 


Pression d'eau interstitielle 


E <=-140 --120 kPa 
o -120 --100 kPa 
o -100 - -80 kPa 
oO -80 - -60 kPa 
O -60 - -40 kPa 
O -40 - -20 kPa 
oO -20 - 0 kPa 

oO 0-20 kPa 

O 20 - 40 kPa 

O 40 - 60 kPa 

O 60 - 60 kPa 

O 8&0 - 100 kPa 
O 100 - 120 kPa 
O 120 - 140 kPa 
O 140 - 160 kPa 


O 160-150 kPa 
E 150 - 200 kPa 
E = 200 kPa 





distance (m) 
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Facteur de sécurité talus amont 


altitude (m) 


o 10 20 30 40 50 60 70 


distance (m) 


Facteur de sécurité talus aval 


altitude (m 


o 10 20 30 40 50 a0 rO 


distance (m) 


20 


20 


100 


100 


AS 


A 
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110 
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b. Cas Retenue vide: 


Facteur de sécurité talus amont 


E 
D 
Le 
5 
T 
distance (m) 
Facteur de sécurité talus aval 
E 
Q 
gë 
3 
6 


distance (m) 


AS 


A 
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altitude (m) 


altitude im) 


c. Cas Retenue remplie à hcritique= 12m : 


La ligne piézométrique 


Pression d'eau interstitielle 


E <-40 --20 KPa 
o -20 - 0 kPa 

o 0-20 kPa 

O 20 - 40 kPa 

O 40 - 60 kPa 

O 60 - 50 kPa 

O 30 - 100 kPa 
oO 100-120 kPa 
O 120 -140 KPa 
140 - 160 KPa 
O 160-180 kPa 
O 180 - 200 kPa 
200 - 220 kPa 
220 - 240 kPa 
O 240 - 260 kPa 
O 260 - 280 KPa 
E 250 - 300 KPa 
E = 300 kPa 


distance (m) 


Facteur de sécurité talus amont 


20 30 40 50 20 vO 


distance (m) 


20 





100 


LE, A 


AS 
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altitude (m) 


altitude (m) 
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Facteur de sécurité talus aval 





30 40 50 La] ro 20 30 100 110 


distance (m) 


d. Cas Retenue remplie à hcritiqgue=5m : 
La ligne piézométrique 


Pression d'eau interstitielle 


m <-100 - -80 kPa 
m -80 - -60 kPa 
o -60 - -40 kPa 
o -40 - -20 kPa 
o -20 - 0 kPa 

o 0-20 kPa 

o 20 - 40 kPa 
o 40 - 60 kPa 
B 60 - 80 kPa 
o 30 - 100 kPa 
O 100 - 120 kPa 
O 120 - 140 kPa 
O 140 - 160 kPa 
oO 160 - 180 kPa 
o 180 - 200 kPa 


200 - 220 kPa 
E 220 - 240 kPa 
E = 240 kPa 





20 30 42 50 60 70 80 30 100 110 


distance (m) 
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altitude (m) 


altitude (m) 





A 
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Facteur de sécurité talus amont 


110 


distance (m) 


Facteur de sécurité talus aval 





distance (m) 
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e. Cas sismique retenue vide : 


Facteur de sécurité talus amont 


altitude (m) 





0 10 20 30 40 50 60 F0 60 90 100 110 


distance (m) 


Facteur de sécurité talus aval 


20 


altitude (m) 


-10 





-20 
0 10 20 30 40 50 60 T0 20 90 100 110 


distance (m) 
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f. Cas sismique retenue hcritique=9 M : 


Facteur de sécurité talus amont 





20 1202 
8 
10 
E 
D 
T 0 
2 
T 
-10 
-20 
0 10 20 30 40 50 60 TO 20 o0 100 110 
distance (m) 
Facteur de sécurité talus aval 
i 1126 
e 


altitude (m) 





0 10 20 30 40 50 60 TÜ 20 90 100 110 


distance (m) 
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g. Les données de SEEP/W pour avoir la ligne piézométrique 


Estimer la fonction de la teneur en eau volumi.. 


Méthode Données de granulométrie 
d'estimation : 


Teneur en eau à 0,5 mžł/m? 


saturation : 
Diamėtre (mm) à 10 % de passant : 


Diamètre (mm) à 60% de passant : 


Limite de liquidité : 
Intervalle de 


SUCOON maximum : 


Nombre de Points : 


Annuler 





Entrer les fonctions de conductivité hydraulique 











Aouter | Supprimer 


conductivité fondation 


Name 1.0e-06 
boue 1.0e-07 
conductivité fondation 1-0e-08 
digue 1.06-09 
drain a 1.0e-10 
È 1011 
Nom: | conductivité fondation Cu 1.0e-12 
= 10613 
Paramètres : Conductivité X vs Pression d'eau interstitielle 5 410e-14 
e 1. 
: p Dre D  10e-1 
Types : Fonction de point de données K hyd i = 
© 1.0e-16fam 
Approx. Exact o 1.0e-17 
Ajustement p 1.0e-18 
de courbe : r 5 K-Saturation : 1 0e-10 
Droit Courbé - Eo 
5 | he 100 To .0e-?0 
1.08-214%— 
Succion capillaire {kPa} Ajouter 
0,01 - ; 
Supprimer 


0,01832980 7 


0,033598183 
0,061568421 
0,11288379 


RS E 


Annuler [+ 





{Fonction de visualisation 
O Afficher la pente de la 
fonction 
(Si Modifier les points de donné 
Ajouter des 
points 


Refaire | Aide 





Succion capillaire (kPa) 





Utiliser l'échelle automatique 


Intervalle de l'axe x : 


Log X 


Log Y 


Plus 


Fermer 





SE, 
LE}, A 
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sal Entrer les Analyses 


Analyses : Ajouter v 


\4 
e A 
LE} | A 


ECOLE MOHAMMADIA D'INGÉNIEURS 








o- geometrie profil dim sud 
5-8 profil dim sud [0 s] 
: [à] Stabilité des pentes [0 s] 
: (à) Stabilité des pentes (2) [0 s] 
EH) profil dim sud (2) [0 s] 
L) Stabilité des pentes (3) [0 s] 
Lf Stabilité des pentes (4) [0 s] 
EH) profil dim sud (3) [0 s] 


[fu] Stabilité des pentes (5) [0 s] 

L Stabilité des pentes (6) [0 s] 

S-E profil dim sud (4) [0 s] 

L Stabilité des pentes (7) [0 s] 

: (à) Stabilité des pentes (8) [0 s] 

EH) profil dim sud (5) [0 s] 

: [là] Stabilité des pentes (10) [0 s] 
: (à) Stabilité des pentes (9) [0 s] 

EH) profil dim sud (6) [0 s] 


fj (Aa) Stabilité des pentes (11) [D s] 








CO X 
Nom : Stabilité des pentes Description : 
Parent : profil dim sud v 
Type d'analyse : Morgenstern-Price v 
Paramètres Surface de glissement Distribution du coef de sécurité Avancé 
Fonction latérale : Fonction semi-sinusoïdale N Valeurs de fn 
Conditions de pression Analyse parente v 
interstitielle de : 


Utilise les résultats de l'analyse parente si c'est SEEP/W, SIGMA/W, QUAKE/W ou VADOSE/W. 


Option d'analyse pseudo-statique par étapes : 


(aucun) v 
Facteurs partiels : 
(aucun) v 








Refaire |v 


Annuler |Y 


Entrer les matériaux 
Materials 


Name 


digue 


drain 


fondation 
géomembrane 


Nom : 








Color 


Affectés… 


Couleur : 


Modèle de matériau : | Saturés / non saturés 








Propriétés 
hydrauliques 
Fn de conductivité 
hyd. : 
Anisotropie — 


conductivité fondation 


Ky'f Ko : 


Vol, Fn de teneur en 
eau ! 

F interst, 
O d'activation : 


fondation 


0 kPa 


soon 









































Annuler |+ Refaire |" 


Afficher la 
légende 


Propriétés... 





130 





LE}, A 


Mémoire du projet de fin d’études =/ AN 


ECOLE MOHAMMADIA D'INGÉNIEURS 


3. Organigramme pour le calcul d'une section rectangulaire d'une poutre en flexion 
simple à l'état limite ultime 


Organigramme pour le calcul d’une section rectangulaire d’une poutre en flexion 
simple à l’état limite ultime 














b „h, fc2s 7 Mu ’ fe 


E 0:85, 


fus : situation accidentelle 


Yb=1.5 : situation durable 


d= 0.9 h 





Oui Non 


Armatures tendues & 
Armatures comprimées 


Mr =u 1 frc.b.d2 
Mr2 = Mr - Mrı 


a ua ES u ı au lieu p u 
M 








s sS Ys 


fe 1.0 : situation accidentelle 





Ys = 1.15 : situation durable 
s: section d’armatures comprimées 
A tion d E 
s: section d’armatures tendues 
As: secti 1 


fe : limite d’élasticité de l'acier 


d'=01h 


Choix des armatures 
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